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Abstract 



We analyze the large-order behaviour in perturbation theory of classes of diagrams with 
an arbitrary number of chains (i.e. photon lines, dressed by vacuum polarization inser- 
tions). We derive explicit formulae for the leading and sub leading divergence as n — oo, 
and a complete result for the vacuum polarization at the next-to- leading order in 1/Nf. 
In general, diagrams with more chains yield stronger divergence. We define an analogue 
of the familiar diagrammatic -R-operation, which extracts ultraviolet renormalon coun- 
terterms as insertions of higher-dimension operators. We then use renormalization group 
equations to sum the leading (lnn/A^/)'^-corrections to all orders in and find the 

asymptotic behaviour in n up to a constant that must be calculated explicitly order by 
order in l/Nj. 
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1 Introduction 



Perturbative calculations in quantum field theories involve integrations over arbitrarily 
small distances and typically lead to divergent results. A finite result is obtained, when 
a theory is first regulatedQ and counterterms are added to the bare Lagrangian Cq. The 
Green functions computed from 

£ = £o(A) + 4?(A) (1-1) 



have perturbative expansions free from divergences order by order in the renormalized 
coupling a{^) so that they are finite in the limit when the cut-off A is removed. A 
justification for this rather ad /ioc-looking procedure can be obtained from the philosophy 
of 'effective field theories': Only the renormalized parameters are accessible to low-energy 
experiments and the apparently divergent regions of integration are in fact (almost) 
insignificant. 

From the beginnings of renormalized quantum field theory it has been recognized that 
the Green functions (in the limit A oo) obtained in this way can not be unambigu- 
ously defined (as certain analytic functions in a neighbourhood of a = 0) through their 
perturbative expansions alone, because they diverge for any a 7^ Although from 
a practical point of view one may consider these expansions as asymptotic (to nature), 
the (non-perturbative) existence of renormalized field theories remains a mathematically 
largely unsolved problem, the divergence of perturbative expansions being one face of 
this problem, the issue of triviality in non-asymptotically free theories being another. 

Without touching the profound problem of existence, the behaviour of perturbative 
expansions as formal series is itself important. In considering the perturbative expansion 
to all orders, one takes in fact a glimpse beyond perturbation theory. Thus, although 
the questions of triviality and Landau poles in general can not be answered without 
knowledge of non-perturbative properties of the theory, some aspects can be investigated 
strictly within perturbation theory. In theories where the coupling can become large at 
low energies, the details of the divergence of the perturbation series may provide some 
hints to selecting the numerically most important corrections already in moderately 
large orders. The dominant source of divergence (at least with present knowledge) was 
identified by Lautrup and 't Hooft , who investigated a particular class of diagrams 
with an arbitrary number of vacuum polarization insertions into a single gauge boson 
line in a loop diagram. The sum of these diagrams has the generic behaviour 



r(p,/x,a(/i)) 



n! n a(/i)' 



n+l 



1 + 01 - 

n 



(1.2) 



Here p is a collection of external momenta and (3o is the first coefficient of the /^-function. 
With the definition given later, the constant a is integer (and sometimes half-integer). 
The ?7,!-behaviour arises as a consequence of renormalization and for this reason has 
become known as renormalon divergence. Note that the term 'divergence' is applied 



-"^Note that there are also schemes without regularization, e.g., BPHZ, differential. Although they 
lead to physically meaningful results there is no bare Lagrangian for them. 
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to the divergence of the perturbation series as well as to the divergences of Feynman 
integrals, which are subtracted in the process of renormalization. 

In strictly renormalizable theories the coupling depends logarithmically on the renor- 
malization scale n and each vacuum polarization loop gives a In k"^/ fi^. This logarithm is 
large whenever the loop momentum k is very different from fi. In the present paper we 
consider only large momentum regions k ^ fi and ultraviolet (UV) renormalons, with 
a > 0. The divergences associated with regions of integration over large loop momenta, 
like d^k/k^ for large fc, are removed by the familiar renormalization procedure. The ultra- 
violet renormalons occur when a large power of In k"^ / //^ is integrated over the remaining 
ultraviolet regions of the subtracted integrands. In particular, the strongest divergence 
(a = 1) is exclusively due to the remaining behaviour d^k/k^ for large k. The fact that 
UV renormalons are related to the large momentum expansion of Feynman integrands is 
crucial for their understanding, since it allows to describe UV renormalon divergence in 
large orders in terms of local operators just as explicit divergences in finite orders in the 
usual framework of renormalization. This observation is the basis for Parisi's hypothesis 
[Q that UV renormalons can be removed by adding higher dimensional operators to the 
Lagrangian. In particular, the leading UV renormalon (a = 1) could be compensated by 
an additional term 

4f = Ae^^(«(^))o- (1.3) 

where the sum runs over all local operators Oi of dimension six. To compensate all UV 
renormalons an infinite series of higher-dimensional operators would have to be added 
to the Lagrangian. In this paper we deal explicitly only with the leading UV renormalon 
and restrict ourselves to operators of dimension six. 

The fact that the removal of ultraviolet renormalons, as the removal of ultraviolet di- 
vergences, can be formulated at the level of counterterms in the Lagrangian implies their 
universality: Once the coefficients Ei have been determined from a suitable set of Green 
functions, the subtraction of the first ultraviolet renormalon is automatic for all Green 
functions. Another consequence is that the Ei satisfy renormalization group equations, 
because Green functions with operator insertions satisfy them. These considerations fix 
the constant h in eq. (|L2|) The solution of the renormalization group equation de- 
pends on a boundary value which remains unconstrained by general considerations and 



is related to the normalization K of the renormalon divergence in eq. ( |1.2| ). 

An alternative approach, based on the Lagrangian at a finite cut-off rather than the 
renormalized Lagrangian, has been taken in The idea is that, since UV renormalons 
arise as a consequence of the infinite cut-off limit, their information is encoded in the 
large-cut-off expansion. Viewed in this way the compensation of UV renormalons bears 
close resemblance to Symanzik improvement of lattice actions . 



Although the absence of the first ultraviolet renormalon as a consequence of eq. ( |1.3| ) 



seems rather obvious from the physical origin of UV renormalons, it has not yet been 
rigorously established. Moreover, the diagrammatic interpretation of eq. ( |1.3[ ) is rather 
unclear: It does not give a clue, which diagrams contribute to the coefficients Ei or 
whether they can be calculated at all in a systematic way. These are the questions which 
we address in this paper. 
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The diagrammatic study of UV renormalons received new attention only recently, 
through the work of Zakharov and others [8-11]. The main difficulty that the di- 
agrammatic approach has to face is that, because the object is to study perturbative 
expansions in large (that is, to all) orders, there is in fact no natural expansion param- 
eter that would select a manageable subset of the infinity of all diagrams. In abelian 
gauge theory it is useful to classify diagrams in terms of complete gauge boson propaga- 
tors. In first approximation, where only the first coefficient in the /3-function is kept, the 
complete photon propagator reduces to a 'chain', a string of fermion loops. With few 
exceptions, previous investigations of UV renormalons have focused on diagrams with a 
single chain. In |[1^ it was shown that at this level one could remove the first ultravio- 
let renormalon from Green functions by counterterms of the form of eq. (|1.3|) . The full 
complexity of calculating the normalization K is already exhibited by diagrams with one 
complete photon propagator: To obtain the value of K, one can not approximate the 
propagator by a string of fermion loops (chain). The exact photon propagator has to be 
kept ||10|, 0. For practical purposes this is equivalent to the statement that K can not 



be calculated exactly. 

An important new insight comes from the work of Vainshtein and Zakharov |l^ 



who investigated the dominant contributions to the large-order behaviour of the photon 
vacuum polarization from diagrams with two chains by making direct use of the fact that 
the UV renormalons originate from the large-momentum regions in loops. The diagrams 
with two chains display a qualitatively new behaviour, because the four-fermion operators 
that appear in eq. ( |1.3D do not contribute to diagrams with a single chain. After insertion 
into the photon vacuum polarization, they were found to yield the dominant large-order 
behaviour. 

In this paper we approach UV renormalons from an entirely diagrammatic perspective 
within the expansion in the number of chains, or, to be precise, in 1/Nf, where Nf is 
the number of fermions. Guided by the interpretation of UV renormalon divergence as 
similar to the usual ultraviolet divergences, we proceed in close analogy with the usual 
renormalization program. We will see that the analysis of UV renormalon divergence 
order by order in the expansion in chains has much in common with the analysis of 
UV divergence order by order in the coupling a. The UV renormalon problem then 
takes a form similar to usual UV divergences: While certain properties like locality of 
counterterms and renormalization group equations can be established to all orders (for 
an arbitrary number of chains), the actual calculation of counterterms is limited to a few 
first terms because of the increasing complexity of integrals. In the present paper we 
proceed heuristically and do not give proofs which could be worked out as generalizations 
of standard cumbersome proofs of renormalization theory. Rather than giving proofs 
we describe the properties of regularized Feynman integrals and subtraction operators 
which would be essential ingredients to these proofs and illustrate how the extraction of 
ultraviolet renormalon divergence works in a number of examples. 

Let us make a remark on the framework of the expansion in the number of chains: 
It might appear inconsequent to replace one expansion (the one in the coupling a) by 
another. Our attitude is that this expansion can shed some light on the organization of 
UV renormalon divergence in the same way as, for instance, the l/A'c-expansion in QCD 
may reveal some information on the strong-coupling regime. In addition, we will see 
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that the dominant UV renormalon divergences in each order of can be resummed 
to all orders in by solving the renormalization group equations. It is possible that 
the conclusions drawn from the l/iVj-expansion are not valid for any finite Nf or valid 
only in a finite range of Nf. However, in an abelian gauge theory we do not consider 
this a likely possibility and rather expect a smooth continuation from the large- iV/ limit 
to small Nf. 

In Sect. 2 we begin with detailing the expansion in chains. We introduce the Borel 
transform as generating function of perturbative coefficients and show how the factorial 
divergence of perturbation theory is encoded in the singularities of analytically regular- 
ized Feynman integrals. We collect some of their properties and classify the subgraphs 
that can contain the dominant UV renormalon. The extraction of the renormalon diver- 
gence then reduces to the extraction of pole parts of analytically regularized Feynman 
integrals at certain positions in regularization parameter space. 

In Sect. 3 we construct an operation that picks out the terms with the largest number 
of singular factors. The leading large-n behaviour at any order in is then found 

by successive extraction of pole parts of one-loop integrals. Beyond the leading large-ra 
behaviour it is essential to apply the method of infrared rearrangement ||15| , p!6| . The op- 
erations introduced in this section are illustrated by the simplest example of the fermion 
self-energy up to two chains. In Sect. 4 we compute the fermion-photon vertex and the 
photon vacuum polarization including all diagrams with two chains and those diagrams 
with three chains that contain an additional fermion loop. 

The universality of UV renormalons becomes most transparent by elevating the di- 
agrammatic subtractions to the level of counterterms in the Lagrangian. We formulate 
the results of the previous sections in this language in Sect. 5. The renormalization group 
equations for Green functions with operator insertions are then used to sum the leading 
UV renormalon singularities to all orders in 1/Nf. We summarize in Sect. 6 and discuss 
further applications of the formalism. 

The derivation of some results quoted in sections 2 and 3 is collected in an appendix. 

2 Renormalons and analytic regularization 

In this section we set up the reorganization of the perturbation series in terms of chains 
and derive the Feynman rules for the Borel transform of this expansion. We show how 
renormalons are related to the singularities of analytically regularized Feynman integrals. 
We consider only the abelian gauge theory (QED) with Lagrangian 

c = --^F^uF^" + v^. + gi^M + 4? • (2.1) 

Since we are interested in large-momentum regions, we can consider the fermions as 
massless. We assume Nf species of fermions, but do not write the flavour index explicitly. 
Summation over repeated flavour indices is understood. The /3-function is given by 

{a = g'/iAn)) 

/3(a)=/x2|^=/3oa' + /3ia3 + ... , /?o = ^ , Pi = ^ ■ (2-2) 
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2.1 Chains 



Consider the perturbative expansion of a truncated Green function G{q;a), where q = 
[qi, . . . iqu) denotes a collection of external momenta and a = a{fi) the renormalized 
coupling. We can organize this expansion in terms of diagrams with complete photon 
propagators 



where n(fc^) is the photon vacuum polarization and ^ the gauge fixing parameter. Each 
such diagram corresponds to a class of diagrams in the usual sense. Let f be such a 
class of diagrams with N complete photon propagators and let the lowest power of g 
that occurs in a diagram in V be g^a^ . (In the abelian theory, e depends only on the 
Green function G and equals the number of external photon lines.) Then we write the 
contribution from F to the series expansion of G as 

oo 
n=N-l 

Note that if we consider a physical quantity, we may take the complete photon propagator 
to be renormalized. Due to the Ward identity, no further renormalization is required 
and the sum over all diagrams with a given number of renormalized photon propagators 
is finite. 

Formally, the series can be written in the Borel representation 

oo 

Gf{q-a) = i-igy^ I dne-"/(^°")i?[Gr](g;n), (2.5) 

where 







oo 

r, 



B[GT]{q:,u)= if/3o-"^" (2-6) 



n=N-l 



is the Borel transform of the series. The factorial divergence of the series then leads to 
singularities of the Borel transform at finite values of u. For example, the large-order 
behaviour of eq. (|1.2|) results in a singularity at m = a and the constant h determines the 
nature of the singularity. The integral in eq. ( |2.5| ) does not exist due to these singularities. 
In the present context, we use the Borel transform only as generating function for the 
perturbative coefficients. We do not consider at all the problem of summation of the 
perturbative series (e.g. by use of the Borel transform). 

Let F be the skeleton diagram corresponding to f . Then Gf is represented as 

„ i4 iV 14, N 

t J i J i J J 
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The momenta pi are assigned to fermion lines and kj to photon hnes. The function /r 
is the Feynman integrand (without the factors of the skeleton diagram, including 

5-functions in momenta from vertices. The complete photon propagators are written 
explicitly, except for the Lorentz structure {—i){g^v — k^ky/k"^), which is included in /p. 
We can drop the piece proportional to ^ in eq. ( ^.31 ) by specifying Landau gauge ^ = 0. 
This will be assumed in the following unless stated otherwise. The convolution theorem 
for Borel transforms can be used to derive the identity 



TV 



B 



n.T 



a 



{u) 



N 



N 



N 



l[duj 



u 



a 



i + n(fc2) 



(u,), (2. 



so that 



B[Gr]{q;u) 



1 



f3S 



N-l 



u 


" N 






llduj 












^ 3 ^ 



N 



]=1 J 



a 



(u,). (2.9) 



This expression is still too complicated, because it contains the complete photon propa- 
gators. The complete propagators can themselves be expanded, if we consider the limit, 
when the number of fermion flavours Nf is large. Then write 



n(fe2 



+ ni(A;2) + — 



Nf 



(2.10) 



This expansion in chains is depicted in Fig. |I|. The dashed line denotes a chain, i.e. a 
photon propagator with an arbitrary number of simple fermion loops inserted. C is a 
subtraction constant for the fermion loop, C = —5/3 in the MS scheme. Then we use 
again the convolution theorem to obtain 



B 



a 



l + n(A;2) 



1 ? ( y 







Hi 



a 



{u-v) + ... (2.11) 



The expansion can be continued to include further terms. We will not need them in this 
paper. The Borel transform of Hi has been computed in [El p[. It can be written as 



B 



Hi 



a 



(u) 
F{u) 



Pi 



, 2 \ " 

-^e^j F{u)-G{u) 



(2.12) 



1 



E 



-irk 



32 



3 n(2 + n) (k^ - (l + u^)'' 



6 



2\ _ 



+ 



+ 2 




+ 



Figure 1: Expansion of the photon vacuum polarization in chains (l/Nf). The dashed 
hne denotes summation over a photon hne with an arbitrary number of fermion loops 
inserted. 



where G{u) is a scheme- dependent function, whose expression, for example in the MS 
scheme, can be found in Sect. 5 of |T^. (Note that the parameter u there, as in most 
publications that deal mainly with infrared renormalons in QCD, is defined with an 
opposite sign compared to the definition we use in this paper.) We will not need an 
explicit form of G{u). 

Let us momentarily replace all complete photon propagators by chains, so that we 
keep only the first term on the right hand side of eq. (|2.11|) . Then 



B[Gr]iq;u) 



N 



Ylduj 



u 



N 

E 



Uj I Gr{q;u) 



(2.13) 



where u = [ui, 



. ,U]\r) and 
GT{q;u) = 



n 



(2vr)4 



N 

n 



(2vr)4 (fe|)i- 



Ir{q;pi,kj 



(2.14) 



is the skeleton diagram with analytically regularized photon propagators.^] Note that 
a different parameter is introduced for each photon line and the fermion lines are not 
regularized. For a physical quantity the sum over all skeleton diagrams is finite. For 
individual diagrams F all UV divergences are explicitly regularized by the parameters m,, 
since the only loops without any regularization parameter in F are fermion loops with 
four or more photon lines attached, which are UV finite. Since fermion lines are not 
analytically regularized, gauge invariance is preserved. 

The factorial divergence of perturbation theory corresponds to singularities of 
B[Gr]{q;u). From eq. ( |2.13| ) together with eq. ( |2.14| ) we deduce that it is sufficient 
to know the singularities of the analytically regularized skeleton diagram. We collect 
results on its singularity structure in the next subsection. 



^Note by comparison with eq. (2.7) that Gr{q',u) contains no powers of the couphng g. 
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Corrections to the approximation of replacing the complete photon propagator by 
a chain can be incorporated through the corrections to the first term in eq. 
B[Gr]{q',u) is then expressed as a certain convolution integral of Gr{q',Ui) and 
B[Ili/a]{u). Similar expressions follow, if one includes yet further terms in eq. 



Thus, provided the vacuum polarization is known to the desired accuracy, one can restrict 
attention to analytically regularized skeleton diagrams. 

When the complete photon propagator is replaced by a chain, the calculation of 
B[Gr]{q',u) for a skeleton diagram with an arbitrary number of chains follows from the 
usual Feynman rules with the following modifications: The photon propagator is given 
by 



S.^--^ -Ta^-' (2-15) 



H) 

(a +2e-prescription is understood) and for internal vertices we haveQ 




/47r 

H)J-^7m- (2.16) 
V Po 



For vertices extending to an external photon line we have 7^ only, cf. eq. ( |2.4| ). Finally 
all Mj-parameters are integrated over with /?o /o" D duj 6{u — J^uj). 

2.2 Singularities of analytically regularized Feynman 
diagrams 

To characterize the singularities of B[Gr]{q;u) we use the results and techniques devel- 



oped for the analysis of divergences and singularities of Feynman integrals [|T8| |I9|, ^0 
Here we summarize the relevant statements regarding their analytical structure. Details 
can be found in the appendix. 

Let FY{q,m) be a Feynman integral corresponding to a graph F. It is a function of 
external momenta q = (gi, . . . , g^) and masses m = (mi, . . . , m^,). Below we really need 
only the pure massless case. The Feynman integral is supposed to be constructed from 
propagators of the form 



Ziip) 



2\l+ri 



(2.17) 



Here Zi is a polynomial of degree ai in the momentum of the l-th line. For the moment we 
let ri be integer. The UV divergences are characterized by the UV degrees of divergence 
of subgraphs 7 of F, 

u;(7) = 4/1(7) - 2L(7) + 0(7) - 2r(7), (2.18) 



■^The factor ^/Att comes from g = V ina, since the Borel transform is taken with respect to a. 
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where h{'y) and L^^y) are respectively the numbers of loops and lines of 7, 0(7) = J^ie-y 

= I]/G7^;- analytically regularized Feynman integral is defined by the replace- 
ment 

+mf ^ + m2)i+n+A,^ (2.19) 

For some lines, the corresponding A-parameters can be equal to zero. 

Let us suppose that there are no infrared (IR) divergences in the graph and that the 
available A-parameters are sufficient to regularize UV divergences in all UV divergent 
subgraphs (in which uj{'y) > 0). The second assumption means that for any such 7 at 
least one of the corresponding regularization parameters is not identically zero. Then 
the analytically regularized Feynman diagram Fr{q;X) can be represented as 

E n (2.20) 

where A = (Ai, . . . , A/,), the functions g^r are analytical in a vicinity of the point A = 0, 
and A(7) = J^ie-y^i- Here the sum is over all maximal forests of F. Remember that a 
forest is a set of non-overlapping subgraphs. A forest is maximal if for any 7 which 
does not belong to JF the set JF U {7} is no longer a forest. The steps that lead to 
eq. ( p. 201 ) are detailed in the appendix. 

We now consider the behaviour of Gr{q',u) of eq. (|2.14|) in the vicinity of the point 
Uo = (""01 ! • • • 5 Moat). The previous discussion is apphcable without modification, when 
all components of Uq are integer. We now allow them to be arbitrary complex numbers. 
Thus we set = —uqi and A; = Uqi — ui. As specific examples, one may have in mind the 
skeleton diagrams shown in Fig. ^ It is natural to define the UV degree of divergence 
of a subgraph 7 dependent on Uq and given by 

^uoi^) = 4/i(7) - Lf{j) - 2Lph{l) + 2no(7) , (2-21) 

where Lf(^ph){''y) is the number of fermion (photon) lines in 7, ^0(7) = J^ie-y^^ (""o/) and 
a similar definition holds for ^(7). Note that eq. ( |2.2(]| ) is derived from the factorization 
of singularities in terms of sector variables, eq. (|A.13|) . The same factorization shows 
that the analytically regularized Feynman diagram Gr{q;u) is a meromorphic function 
of u with poles described by the equations 

n(7) = no(7)-A(7) = -M7)/2]+fc, A: = 0, 1, 2, . . . , (2.22) 

where [x] denotes the integer part of x. Since A (7) is by definition small, we see that 
the subgraph 7 contributes a singular factor 1/A(7) to eq. (|2.2CI| ) only if ^0(7) is integer. 
Thus Gr{q',u) can be represented as 
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Figure 2: Ladder contributions to the photon vacuum polarization. 



Gr{q;u)=J2 n —r\ ^-^9A2^) (2.23) 

V . ^ \ ^>n ^0(7) -^i(7) 
1^0(7) integer 

where the functions gjr are analytical in a vicinity of the point Mq. Note that because 
the fermion loop with four photon lines attached is UV convergent in QED, the forests 
that contain this subgraph as their smallest element should be excepted from the sum 
above. 

When the singularities of Gr{q;u) at ^(7) = are removed by the usual process of 
renormalization, the condition that ^(7) is integer for at least one 7 G in order to 
obtain a singularity implies that the Borel transform B[Gr]{q', u) is analytic for < m < 1. 
Moreover, the singularity closest to the origin at m = 1 arises only from the boundaries 
of the integration over the Uj in eq. (|2.13|) . The singularity at m = 1 is called the 



leading ultraviolet renormalon. The strongest singularity for a given set of diagrams 
r comes from those forests which contain the maximal number of divergent subgraphs 
{uju^{'~f) > 0) and from those points Uq in the integration domain, where ^0(7) is integer 
for all divergent subgraphs 7 of the forest. 

As long as we are interested only in the leading UV renormalon, we are interested 
only in the domain < Ui < 1, < J^ui < I, whence ^0(7) < 1- Therefore ujug{'j) > 
implies uj{'y) > —2 (here 0^(7) denotes the usual degree of divergence, i.e. ^0(7) = 
in eq. ( p.21|) ). In particular, if a subgraph 7 includes all the regularized photon lines, 
we have exactly ^^u^i'^) = ^^il) — 2. Thus we have to consider not only divergent (in 
the usual sense, i.e. with u^j) > 0) subgraphs but also those convergent ones, that 
have Lv{'y) = —1 or (^(7) = —2. An example of such a subgraph is given by the 'box', 
consisting of lines {1, 2, 4, 7} in Fig. |^. Thus the IPI Green functions with the following 
field content contain the leading UV renormalon: ipip, tpipiljilj , ipA^ip {n = 1,2,3). All 
other Green functions are superficially free from the leading UV renormalon and develop 
it only through subgraphs that can be classified in terms of the Green functions listed 
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above. 

Barring cancellations, eq. (|2.23|) allows us to derive in a straightforward manner the 
n-dependence of the large-order behaviour of a set of diagrams T. As an illustration, 
we consider the contributions to the photon vacuum polarization depicted in Fig. ^. 
Referring to diagram (b), denote 71 = {1, 3, 6}, 72 = {1, 2, 4, 7}, 7+ the union of 71 and 
72 and r the entire graph. Thus u(7i) = ui and ^(72) = ^(74.) = u(T) = ui + U2 = u. 
The last equality follows from the delta-function in eq. ( |2.13| ). Let JF^ = {71, 7+, F} and 
JF2 = {72,7+,F}. When u approaches unity, we find that in the vicinity of the point 
'1,0) the singular factors in eq. ( |2.23D for these two forests are given by 



Tx: V^i ^ -^2: TT-^- (2.24) 

(1 — uy \ — u\ (1 — uY 

The integration over u\ and U2 is trivial and leads to the singularities of B\Gy\{(i\u) 



Tx: ^"^^ ^ i^i/J^lnlnno'^+i (2.25) 

(1 — uy 



T2: ^ i^2/?o^!n2a-+i. (2.26) 

(1 — u)-^ 

The constants Kx 2 will be determined in Sect. 4. To the right of the arrows we 
have indicated the corresponding large-order contribution to the perturbation expan- 
sion {—igY'Yli'^n(y^^^ of the vacuum polarization, which can be deduced from eq. ( p.6|) . 

We notice that the strongest divergence, n\r? ^ arises from the forest that contains 
the box subgraph. The contribution from one chain, shown in Fig. |]a, produces only 
1/(1 — n)^ P]. Therefore the diagrams with two chains dominate by a factor of n over 



the single chain ones [jT4|. This enhancement has a simple interpretation in terms of 
counting logarithms that are integrated over. In a given order g^a^^"^ in perturbation 
theory the single chain in diagram (a) gives n logarithms from n fermion loops. One 
additional logarithm is generated from the UV behaviour of the reduced diagram, when 
the subgraph {1, 2, 4} is contracted. The result is n! n. In case of (b), at the same order 
in perturbation theory, we have only n — 1 fermion loops, but the reduced diagram after 
contraction of the box subgraph contains two UV logarithms so that the total number 
of logarithms is the same as for (a). However, there exist n ways of distributing n — 1 
fermion loops over two photon propagators and the stronger divergence, n! arises from 
this combinatorial enhancement. In the following we will see that the singular factors 
in eq. ( p.23[ ) originate from d'^k/k^ terms in the expansion of Feynman integrands in 
external (exceptional) momenta. This will allow us to associate the box subgraph with a 
counterterm proportional to a four-fermion operator in the sum of eq. ( |1.3|) and to give 
yet another interpretation of this enhancement. 

What can be expected from diagrams with more than two chains? Consider the 
ladder diagram with chains in Fig. It is not difficult to see that the strongest 
singularity at m = 1 comes from a forest built as follows: Start with the subgraph 
{1,2, + 2, 2A^ + 3} and continue by including subsequent ladders to the right. The 
product of A^ -|- 1 singular factors for Mq = (1, 0, . . . , 0) is given by 
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resulting in the singularity ln^~^(l — «)/(! — u)^. Thus each additional chain yields an 
enhancement only by a factor of Inn for large n, but no additional factors n occur. If we 
anticipate the association of the box subgraph with a four-fermion operator insertion, 
we can interpret the dressing by ladders as a renormalization of this operator. Thus we 
expect that the series of singularities generated by ladder diagrams (and many others) 
can be summed by renormalization group methods and exponentiates according to 

This will be discussed further in Sect. 5. 

Beyond the point u = 1, the Borel transform B[Gr]{q]u) defines a multi- valued 
function due to the cuts attached to the singular points. The above analytic properties 
of Gr{q', m) imply that after integration over u the singular points in u in the right half- 
plane occur at integers with a cut attached to each such point. We can therefore conclude 
that to any finite order in the expansion in chains the only singular points in the right half 
of the Borel plane are UV renormalons at integer u. Whether this reflects the correct 
singularity structure of QED, depends on the behaviour of the expansion. For 

example, the number of skeleton diagrams grows rapidly in higher orders of 1/iV/, so that 
the 1/A^j-expansion of the normalization K of UV renormalons could be combinatorically 
divergent. Factorial divergence of the perturbative expansion due to the number of 
diagrams is indeed expected for theories with bosonic self-interaction. For theories with 
no bosonic self-interaction such as QED, the Pauli exclusion principle enforces strong 
cancellations. For QED with finite UV cut-off, the combinatorial divergence was found 



to be 21, 22 



rl^j (-)„'. (2.29) 



If this combinatorial behaviour persists when it interferes with UV renormalization, the 
corresponding singularities occur at |m| = oo in the Borel plane. Thus it is reasonable to 
assume that in QED the conclusions obtained from the chain expansion pertain to the 
full theory in any finite domain in the Borel plane. 

It is interesting to compare this with the non-abelian gauge theory. In this case 
the existence of instantons and bosonic self-interaction leads to singularities at finite 
values of u, connected with the value of the action of an instanton-antiinstanton pair 
[ p3| . Because the action is independent of Nf, this singularity does not show up in 
any finite order in the 1/iV/ expansion. (The simplest way to see this is to rescale the 
coupling as a = aNf and to write exp{—S/a) = exp{—NfS/a), which has vanishing 
Taylor expansion in 1/Nf.) Thus, as physically expected (for other reasons as well), a 
strict expansion is certainly in trouble for non-abelian theories. 
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3 Extraction of singularities 



While the formulas given in the previous section allow us to derive the n- dependence 
of the large-order behaviour for a given set of diagrams, we are still lacking a method 
to compute the overall normalization without having to compute the analytically regu- 
larized Feynman integrals exactly. In this section we first derive a formula that allows 
calculation of the leading singularity, defined as the collection of terms with the maximal 
number of singular factors in eq. ( |2.23| ), by consecutive extraction of pole parts of one- 
loop integrals. The next-to-leading singularity is defined by the collection of terms with 
one singular factor less than the maximal number. We shall see that the computation 
of the first correction to the normalization K of the large-order behaviour requires some 
specific contributions to the next-to-leading singularity and corrections to the approxi- 
mation of the complete photon propagator by a chain. The techniques developed in this 
section are applied to the fermion self-energy for illustration. 

3.1 The basic formula for the leading singularity 

As mentioned above the leading singularity (LS) of the given Feynman integral is defined 
as the sum of terms in eq. (|2.2CI|) with the maximal number of the factors A (7) in the 
denominator. Eq. ( |2.23| ) shows that this number is equal to or less than the maximal 
number of divergent subgraphs that can belong to the same forest. Note that 'divergent' 
for given Uq means uju^i^j) > 0. Suppose that these divergent subgraphs form a nested 
sequence 71 C . . . C ■jno and let ki, i = 1, . . . ,no, be the loop momentum of 7j/7i_i 
(70 is defined to be the empty set). The singular factors in eq. ( |2.23| ) arise when the 
internal momenta of a given subgraph are much larger than the external momenta of 
this subgraph. Therefore we expect that the leading singularity arises from the strongly 
ordered region 

/ci > A;2 > . . . > A;„o . (3.1) 

Because of this ordering every subgraph appears as insertion of a local vertex with respect 
to the next loop integration. The leading singularity can then be found from consecutive 
contraction of one-loop subgraphs and insertion of the polynomial in external momenta 
associated with the local vertex. This fact is succinctly expressed by the following simple 
representation for the leading singularity: 

LS{Fr)= n ^A^^^^x^P-yh-i^^))) ■ (3-2) 

^■\^hMv\=no'y&^h,div 

Recall that A (7) = ^0(7) — ^(7) and that for the given point Uq the values ^0(7) must 
be integer for all 7 G J^h,div to obtain the maximal number of singular factors. Here each 
maximal forest is represented as J-'hU J-'r = J-'h,div U J-'h,conv U where the subscripts 
'div' and 'conv' denote IPI elements respectively with uju^{'j) > and u^ugil) < 
and J-'r contains all non-lPI elements. The number no denotes maxjr i.e. the 

maximal possible number of divergent subgraphs that can belong to the same maximal 
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forest. Furthermore 7_ is the set of maximal elements 7' G with 7' C 7. Each factor 
resA(7)-F-y/^_ (m(7)) is a polynomial with respect to external momenta and internal masses 
of 7/7. . It is implied that these factors are partially ordered and before calculation 
of the residue resA(7)-F7/7_ (^(7)) all polynomials associated with the set 7_ are inserted 
into this 'next' reduced diagram. Note also that for the Feynman integral F^/7_(u(7)) 
the sum of regularization parameters is m(7), rather than u(7/7_). It does not matter 
into which line of the reduced diagram 7/7. the regularization parameter A = ^(7) is 
introduced, since the corresponding pole part in this A does not depend on this choice. 
Eq. ( ^.2| ) formalizes the expectation that the leading singularities of the Borel transform 
can be calculated by extracting pole parts of one-loop subgraphs and reduces to the 
purely combinatorial problem of writing down all maximal forests for a given diagram 

r. 

A proof of eq. (p.2|) , based on the a-representation technique to resolve the singular- 
ities of Feynman diagrams, can be found in the appendix. Here we summarize only the 
main points of this approach. 

The initial step is to represent the Feynman diagram as an integral over L positive 
parameters ai corresponding to its lines. Then one performs a decomposition of this a- 
representation into subdomains which are called sectoriQ and correspond directly to one- 
particle-irreducible subgraphs of the given graph. After introducing, in each sector, new 
variables associated with the family of IPI subgraphs of the given sector, the complicated 
structure of the integrand is greatly simplified, and the analysis of convergence and/or 
analytical structure with respect to the parameters of analytic regularization reduces 
to power counting in one-dimensional integrals over sector variables. At this point one 
observes that the singular factors exactly correspond to divergent subgraphs of the given 
graph. 

The leading singularity then appears from the sectors with maximal number of di- 
vergent subgraphs. To calculate coefficients of the products of these singular factors 
one uses the local nature of UV divergences. Practically, this means that calculation of 
residues with respect to the corresponding linear combination of analytical parameters 
reduces to Taylor expansion of the corresponding subdiagram in its external momenta 
and inserting the resulting polynomial into the reduced diagram. Eventually one comes 
to eq. Q. 

An alternative proof could use0 the method of glueing [^] which rests on integra- 
tion of a given diagram with an additional (glueing) analytically regularized propagator. 
Within this method, the information about the large momentum behaviour is encoded 
in analytical properties of glued diagrams with respect to parameters of analytical reg- 
ularization. In our problem, the idea would be to use an inverse translation of these 
properties to get the analytical properties of a diagram in the u-parameters from the 
large momentum expansion of specific subdiagrams. 



^ Eq. (3.1) provides an example ol such a sector in the momentum space language. While this 
language is useful for heuristic arguments, the a-parametric technique is adequate for proofs because of 
the simplicity of singularities in terms of sectors and sector variables. 

^We are grateful to K.G. Chetyrkin pointing out the possibility of such a strategy. 
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Figure 3: Self-energy diagrams at next-to-leading order in 1/Nf ((c)-(e)). (a) and (b) 
show the non- vanishing subgraphs for diagrams (c), (d). (f) and (g) show the reduced 
diagrams with insertion of a polynomial in external momenta of the contracted subgraphs 
'c) and (d). 
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3.2 Self-energy: The leading singularity 

In this subsection we exemplify the extraction of the leading singular behaviour by self- 
energy diagrams with two chains. These diagrams are shown in Fig. ^ and d. Diagram 
e also contributes at the same order in 1/Nf. (For the counting in Nf it is useful to think 
of a as being of order 1/Nf.) It corresponds to a correction to the photon propagator 
in eq. ( |2.11| ) and will be dealt with later. For completeness, we note that the Borel 
transform of the leading order diagram Fig. ^ is given by 




(3.3) 



The self-energy requires renormalization beyond renormalization of the coupling. In 
general this requirement is translated into the necessity to subtract those singularities in 
the Uj that give rise to a singularity at m = of the Borel transform. (Such a singularity 
prevents the calculation of coefficients r„ of the perturbative expansion as derivatives 
of the Borel transform at m = 0, see Sect. 2.1.) In the present case, the pole at u = 
is absent, because the logarithmic UV divergence of the one-loop self-energy cancels in 
Landau gauge. When fermion loops are inserted in the photon line, logarithmic overall 
divergences are present. In terms of the Borel transform their subtraction in a specified 
scheme amounts to subtracting an arbitrary function of u with the only restriction that 
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it does not have a pole at m = 0. In the MS scheme this function is entire |11] and does 
not introduce new singularities at any u. In general, it is quite non-trivial to determine 
this function in a given renormalization scheme (for a single chain see [|11|], appendix A). 
We do not discuss this point further, since our main interest concerns physical quantities. 
In this case all subtractions necessary for UV finiteness are implicitly contained in the 
definition of the renormalized complete photon propagator. 

Let us now turn to the diagrams c and d in Fig. ^. Each maximal forest contains two 
non-trivial elements, the entire graph F and a one-loop subgraph 7. The contributions 
from most forests vanish, because the reduced graph F/7 is a tadpole graph, for instance 
in case of 7 = {1,2,4,5} for diagram c. For diagram c only one non- vanishing forest 
remains, J-' = {{2, 4}, Fc}, for diagram d we have {{1, 3, 4}, F^} and the symmetric forest 
{{2, 4, 5}, F^}. The leading singularity has two singular factors: 1/(1 — m) from F/7 and 
1/(1 — Mj) (j = 1,2) or 1/uj from 7. Thus the leading singularity at m = 1 arises from 
the points Uq = (1, 0) and Uq = (0, 1) in regularization parameter space. 

We consider first diagram c and the vicinity of (0, 1). If we denote the momentum 
of lines 2 and 5 hj Q — ki, the pole part of 7 = {2,4} is given by the second line of 
eq. ( |3.3| ) with u replaced by U2 and Q hj Q — ki. The insertion of the polynomial in 
Q — ki associated with the contraction of 7 into the reduced graph is shown in Fig. 
The corresponding Feynman integral (up to constants that can be read off eq. (|3.3|) ) is 
given by 



1 r d%j^{fQ-h){Q-h?{^-hm-hh. 1 { k^K\ 
I-U2J {2n)^ {Q-k,f {klY-^V'"'' J ■ ^■'■^^ 

The denominator in Q — ki is cancelled and the integral vanishes. In the vicinity of 
the other point, Ui=l, U2 = 0, the potential pole l/u2 is absent due to finiteness of 
the one-loop self-energy in Landau gauge as noted above. Thus diagram c does not 
contribute at all to the leading singularity. Note that according to eq. ( p.2|) and the 
appendix, the regularization parameter associated with the photon line of the reduced 
diagram is u = ui + U2, rather than ui. In other words, the reduced diagram inherits 
the M-parameters from the contracted subgraph, here U2- The reason for this is that 
the singularities in u from an element 7+ G (in the present example 7+ = F) are 
determined by ^(7+), rather than m(7_|./7), see Sect. 2.2. A less formal argument comes 
from doing the subgraph 7 exactly. Simply by dimensional reasons the Borel transform 
is proportional to (yU^/ {Q — ki)"^)'^^, so that the parameter U2 is attached to the fermion 
line of the reduced graph. The singularities of the reduced graph arise in turn from the 
integration region ki ^ Q, so that we may neglect Q in this expression and attach the 
parameter U2 to the photon line, where it combines with ui to u. This statement is true 
in general (see the appendix): In calculating the leading singularity, the w-parameters 
inherited from the contracted subgraph can be attached to an arbitrary line of the 
reduced graph. 

Turning to diagram d, consider the forest {{1, 3, 4}, F^}. Again we have to analyze 
the points Mo = (1; 0) and Uq = (0, 1). The second one does not contribute two singular 
factors, since the one-loop vertex subgraph {1,3,4} shown in Fig. |^ is UV finite in 
Landau gauge, so that no singularity 1/ui occurs. In the vicinity of the first point 
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we require the singularity of the vertex subgraph at Mi = 1. With the momentum 
assignments of Fig. ^ we find 



(2^)4(p-A;)2(p'-A;)2(A;2)i-«i 
-.1 i 1 



167r2 1 — ui 



)4 



(3.5) 



where Q = p' — p is the momentum of the external photon. Note that the terms not 
containing Q in the second line vanish on-shell and in the Lorentz structure of the terms 
containing Q one recognizes the Feynman rule for insertion of the operator ip'-jyipd^^F^^^ . 

In order to calculate the pole parts of one-loop integrals like the vertex above, one 
should use the fact that the poles of interest are of ultraviolet origin, that is from the 
region of integration, where k is much larger than all external momenta. Therefore we 
can expand the denominator in p/k and p' /k. The IR divergences that arise in this 
way can be regularized by a cut-off \k\ > /i or in any other convenient way. Lorentz 
invariance then allows us to drop all terms with an odd number of /c's in the numerator 
and to simplify the numerator by relations such as kpk^ — > Qpa/^ and its generalizations 
to more factors of A;. If the original integral was logarithmically UV divergent, the result 
is expressed as the sum 



n=0 nl.. 



|fc|>/i 



(2^)4 (fc2)2+n-«(7) ' 



(3.6) 



where Vn{<i) is a polynomial of degree n in the external momenta. The A;- integral gives 
a pure pole, so that Vn{c[) can be identified with the residue in eq. ( p.2|) that is to be 
inserted into the reduced graph. For the fermion-photon vertex Vo{q) vanishes and for 
Pi(g) we obtain eq. Notice that the pole at m = n is related to the term d^k/ (/c^)^"*"" 

in the expansion of the integrand. The relation to Taylor operators is clarified further 
in the appendix. 

The remainder of the calculation is straightforward. Eq. (|3.5|) is inserted into the 
reduced graph shown in Fig. |^ and the singular part of the reduced graph is extracted 
in the same way as before. Adding the contribution from the symmetric forest, we arrive 
at 



G'rd(Q;'"i,'"2) 



(4vr) 



1 



1 — Ul — U2 



1 — Ul 



+ 



U2 



(3.7) 



After integration over ui and U2 according to eq. (|2.13|) , the final result for the leading 
singularity aX u = 1 from the two-chain diagrams c and d is 



SrJS](n)+i?r,P](7 



/.2 



1 1 



2 1 - u iV/ 



[In(l - n) + 0{\)\ 



(3.8) 
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which has to be compared with the leading order contribution, eq. ( p.3| ). (The factor 
originates from 1/ (3q in eq. ( p.l3|) .) In this case we find (S^nl Inn for large n with 
a logarithmic enhancement as compared to the leading order. The 0{1) contribution in 
brackets will be found in Sect. 3.4. This concludes the illustration of the general method 
of extracting the leading singularity. The algebraically more extensive case of the vertex 
and vacuum polarization is treated in Sect. 4. 



3.3 Next-to-leading singularity and IR rearrangement 

At this point it is helpful to refer again to the large-order behaviour of perturbative 
expansions, eq. ( |1.2| ). Within the expansion in the number of chains, we write 

K = K^'^\l + 5K^'^^ + ...), 6 = 6™ + ... , (3.9) 



where 5K^^\ 6'^' denote contributions from N chains.^ Since we expand in chains before 
taking the large-n limit, eq. ( p..2| ) must be written as 



r„ = [1] n! n^''' (l + {6[2l In n + 6K^^^ } + ...) 



i + o[ - 

n 



(3.10) 



where the ellipses represent terms like 6^^]^ In^ n etc. from diagrams with three and more 
chains. Close to = 1, the corresponding Borel transform eq. ( p.6|) is 



B[G]{q;u) 



i^Mr(i + 6M) 
(1 _ uY+hW 

X "l + C'((l-'u)ln'=(l-n) 



(l + {-6[2l (ln(l-u)-V'(l + 6W)) +<5K[21|... 



(3.11) 



where il){x) is Euler's T/^-f unction. From this expression we deduce that 6'^' can be read 
off from the leading singularities of diagrams with two chains, such as those considered 
in the previous subsection. On the other hand to obtain the first correction SK^"^^ to 
the normalization, we require also some terms with only one singular factor. This can 
be either 1/(1 — ui — U2) or 1/(1 — Uj) {uj = 1,2). After integration over ui, U2, the 
second type produces a singularity ln(l — u) and contributes only to 1/n-corrections in 
eq. ( ^.1(J| ). Thus, we consider only singular factors of type 1/(1 — -u). 



This distinction is important both as far as the physical origin of the singularity is 
concerned as well as its extraction. Recall that, for example for the contribution from 
Fig. ^ to the self-energy, the leading singularity (two singular factors) comes from the 
loop momentum regions 

Q<^ki<^k2 Q<^k2<^ki. (3.12) 

To obtain at least one singular factor, we must consider the regions 

^To be precise, different K and b should be introduced for each operator that appears in eq. (pT^). 
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Q <. ki ^ k2 , 

Q ~ /ci < /C2 Q^k2<^ki 



(3.13) 



Only in the first case we obtain 1/(1 — m), because both virtual photons have momentum 
larger than the external momentum Q, so that the large-momentum part of the diagram 
contains all available w-parameters. Moreover, because of this, the residue of the pole 
is polynomial in the external momentum. On the other hand, the regions listed in the 
second line lead to singularities 1/(1 — Uj). The residue is no longer 'local' (polynomial 
in external momenta), but contains ln(Q^//i^). A simple way to see this is to note that 
by dimensional reasons the Borel transform for diagram d is proportional to Q (Q^/ /i^)"- 
This implies that the leading singularity at m = 1 in fact arises in the combination 



1 , Q' 

In ^ 

— u /x^ 



(3.14) 



so that the coefficient of the logarithm is related to the residue for the leading singularity. 
One may compare this with the 1/e^ and (1/e) ln(Q^//i^) poles for a dimensionally regu- 
larized two-loop integral before the subtraction of sub divergences is taken into account. 

We also see that if explicit renormalization (in the usual sense) is needed as in case of 
the self-energy, the arbitrariness in choosing finite subtractions affects only the residues 
of 1/(1 — Uj). Indeed, the effect of such a subtraction is described by (taking j = 2, so 
that the logarithmic UV divergence occurs in the subgraph with photon line regularized 
by Ml) 



1 



1 — UUl 1 — U2 



- + G(ni) 

Ml 



(l-M)(l-n2) 1 - n2 



°<'"' (3.15) 



where G{ui) has no pole at ui = 0, but is arbitrary otherwise. Recall that a is in fact 
zero for the self-energy (in Landau gauge). We therefore conclude that the terms that 
contribute to SK^"^^ are local and do not depend on the renormalization scheme. 

The fact that the residues of the singular factors that contribute to the overall nor- 
malization of renormalon divergence are polynomial in the external momenta allows us 
to use the method of IR rearrangement [|15], [1^ (see also [^] for a review) which was 



originally developed for and successfully applied to renormalization group calculations. 
In its original form, the method is based on polynomial dependence of UV counterterms 
on momenta and masses. It consists of the following two steps: (a) Differentiation in 
external momenta and internal masses. The problem thereby reduces to the calculation 
of zero-degree polynomials, (b) Putting to zero all the internal masses and external 
momenta except one. Usually, at step b one puts to zero all the momenta and masses 
and chooses a new momentum that flows through the diagram in some appropriate way. 
The goal is to make this choice in such a way that the diagram becomes calculable. In 
simple cases the problem reduces to calculation of the pole part of primitively divergent 
propagator-type integrals or vacuum integrals with one non-zero mass. The application 
in the present context is illustrated in Sect. 3.4. 
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The arguments presented above can easily be generalized to the contribution 5K^^^ 
to the normalization K from chains. For such a contribution we need at least one[| 
singular factor of type 1/(1 — m), (with u = Ui + . . . + Un). Although no ordering of 
loop momenta is required, all loop momenta ki have to be much larger than the external 
momenta q in order to get such a singular factor. Therefore the residue of these terms 
remains local for any A^. Similarly, if UV renormalization is necessary, the arbitrariness 
in choosing finite subtractions affects the singularity of the Borel transform only at 
the level of ln'^(l — u) (where typically k < N) and therefore modifies only the 1/n- 
corrections to the asymptotic behaviour in eq. ( |3.10| ). A scheme-dependence enters the 
overall normalization only through the subtraction constant C for the fermion loop p6| 



see eq. ( p. 10 ). With the help of IR rearrangement, the calculation of K to an arbitrary 
order in the expansion in chains reduces to extracting pole parts of primitively divergent 
propagator-type integrals or vacuum integrals with one non-zero mass. 

Note that IR rearrangement can be applied also to the calculation of 1/n corrections 
to the asymptotic behaviour. Because of the logarithms in external momenta that enter 
in this place, one must apply explicit subtractions in subgraphs, which remove these 
logarithms. Since the combinatorial structure of eq. (|2.23| ) is completely analogous to 
the one that arises in the context of usual renormalization, the combinatorial structure 
of the subtraction operator is equally simple. We will return to this point in Sect. 5. 



3.4 Self-energy: The next-to-leading singularity 

In this subsection we illustrate the calculation of 6K^'^^ and the use of IR rearrangement 
for this purpose by continuing with the fermion self-energy. When F = Fc, F^ as shown 
in Fig. |c and d, we have for the integral in eq. ( p.l3|) the following asymptotic behaviour 
at M — > 1: 



r du' GriQ; u, u - u') [ai{Q) ln(l - u) + a2{Q) + O {{I - u) ln(l - u))] , (3.16) 

JO I — u 

where are proportional to QQ^ with no logarithm in j [J^ as explained in Sect. 3.3. 
(ai(Q) is given in eq. ( p.8|) .) According to eq. ( p.23|) , the function Gy{u\,U2) can be 
represented as 



Gr(ui,U2) 



9l{Ul,U2) _^ g2{Ul,U2) 
1 — til 1 — U2 



h{ui,U2) 



(1 - til - U2){1 - Ui){l - U2) 



(3.17) 



Now we apply the following proposition: If f{u') is analytic at < m' < 1 then 



du' - — 

/o {I - u'){l - u + u') 
when u ^ 1. As a result we have 



[/(0) + /(l)]ln(l-tx) 



(3.18) 



""if fel^l is different from zero (as for the vacuum polarization, see Sect. 4), 6'^' + 1 such factors are 
needed. 
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h{ui,U2) 



rU 

/ dui<lu2 6{ui+U2-u)- 

Jo [l - Ui - U2){1 - Ui){l - U2) 

= 0) + /i(0, 1)] + r du'hiu', 1 - u'), (3.19) 

1 — n 1 — ti Jo 



when u ^ 1. Here 

- 1 - u') h{l,0) h{0,l) 

h(u , 1 — n ) = ; ; 

^ ' ' u'{l-u') l-u' u' 

giiu',l-u')-giil,0) ^ g2iu',l-u')-g2{0,l) ^o) 
1 — u' u' 

Define bi = gi{l, 0) and 62 = fi'2(0, 1), so that the leading singularity is determined by 
Gr{Q;ui,U2) = — near ni = l,n2 = 0, 

L — Ul — U2 1 — Ui 
1 62 

Gr{Q;ui,U2) = near ui=0,U2 = l. (3-21) 

1 — Ul — U2 1 — U2 

Then the coefficients that enter eq. ( p.l6| ) are expressed as 
ai = -bi - 62, 



a2 = [ du' \veSu-^iGr{Q;u',u- u') - -^\ 
Jo I I -u' u' } 



(3.22) 



Generally speaking, one can directly use eq. ( p.22|) for calculation. This is indeed 



possible for the rainbow diagram of Fig. |^c, which can be computed exactly for arbitrary 
Ul and U2. After a straightforward calculation, we get 



ai 

h{u , 1 — u') 



(see Sect. 3.2) 



(4vr) 



3\ l-u' 



2 J 2 + u' 



(3.23) 



Integration according to eq. (p.22|) and reinstating the factors from eq. ( p.l3|) yields 



in 



1 1 



1 



21-u N 



f 



27 2 9 
— In - + - 
4 3 4 



(3.24) 



In most cases, however, the analytically regularized Feynman integral is hardly cal- 
culable for arbitrary values of ui and U2- This is true in particular already for diagram d 
in Fig. ^. In this situation one ffist calculates coefficients bi and 62 (and if necessary also 
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the coefficients of terms 1/((1 — u)uj)) with the help of general statements regarding the 
leading singularity, and then the function that enters the integrand of eq. ( p.22|) using 
IR rearrangement. To apply IR rearrangement to diagram d we use the fact that the 
coefficients are polynomials in Q, in this case QQ^. Thus we differentiate three times. 
Since Oq QQ"^ = 48, the substitution 



Gr{Q;u',u 



m^^dl^^QGr{Q;u,u-u') 



(3.25) 



leaves the coefficient 02 (Q) unchanged. To reduce the number of terms that arise in the 
process of differentiation, we route the external momentum Q through the lines 1 and 5 
in Fig. I^d, so that the corresponding momenta are Q — ki and Q — k2, respectively. The 
threefold differentiation of the product of photon propagator and fermion propagator 
gives rise to six new diagrams Fj with the same topology as the original diagram but 
with propagators differentiated in a special way. The residue 



ieSu-iGriiQ;u' ,u 



(3.26) 



for each diagram is now a number independent of the external momentum Q. The second 
step in the IR rearrangement is to apply the equation 



res„_iGr, {Q;u,u-u)= res„_iGr^ {p; u, u 



u 



(3.27) 



The meaning of this equation is as follows: Instead of the initial external momentum Q, 
one introduces a new external momentum p which flows through the diagram in some 
appropriate way. In our example diagram d (and its descendents Fj) we choose p to flow 
through line 5 (so that the vertex to which lines 1, 4, 5 are attached becomes external). 
Due to this choice of external momentum the 'new' diagram G'-p.{p;u' ,u — u') becomes 
explicitly calculable by consecutive application of one-loop integrations. 

The extensive algebra connected with numerators of integrals after differentiation 



can be done with the help of REDUCE A generic integral has the form 



N{ki,k2,p) 



(27r)4 (27r)2 {klY{{ki - k^f fiklYHp - ^2)^) 



(3.28) 



We can take the /ci-integral. Since the remaining k2 satisfies k2^ p to obtain a singu- 
larity 1/(1 — m), we can then extract the pole factor from the fc2-integral by expansion 
in the external momentum as in Sect. 3.2. Summing over all six contributions from 
differentiation, we obtain 



ai 

h{u\ 1 — u') 



(4vr) 



(see Sect. 3.2) 

1 1 1 
+ 7 + 



1 



2 + u' 1 + n' 2-u' 3-n' 



(3.29) 



Integration according to eq. (p.22| ) and reinstating the factors from eq. ( p.l3| ) yields 
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1 1 



1 



21-n 



/ 



ln(l - u) + 21n2 - ln3 + - 



(3.30) 



When combined with eq. (|3.24|), this results in 



5Kfl^ = j(351n2-311n3 + 12) . 



(3.31) 



3.5 Self-energy: Propagator corrections and summary 

To complete the calculation of all contributions at subleading order in 1/Nf (for the 
A^j-counting one rescales a = aNf) we have to analyze diagram e in Fig. ^ (plus the 
self-energy type contributions to the vacuum polarization insertion, not shown in the 
figure). This diagram is treated differently from c and d, because it appears as the first 
correction to approximating the complete photon propagator in eq. ( p.ll| ) by a chain, 
the first term on the right hand side of eq. (|2.11| ). Corrections of this type have been 
considered previously in To include the second term in eq. ( |2.11| ), we write, as 

in 
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reg 



F{u) + 



u 



Grcg(n) = G{u) + 



(3.32) 



where F{u) is given in eq. (|2.12|) and Fj.eg{u) and Greg(w) are finite at m = 0. When 
the second term on the right hand side of eq. ( |2.11| ) is inserted for the single complete 
photon propagator in eq. (|2.9|), we obtain 



Po 



^ J dvv{B^^[mu)F,,g{u-v)-B^jmv)G,,g{u-v)) , (3.33) 



Pi 



dv V 



u — V 



where i?f^[E](M) is the leading order contribution given in eq. ( |3.3| ). The vacuum polar- 
ization insertion requires renormalization and the function G{u) takes into account the 
arbitrariness of finite subtractions, that is, the arbitrariness in defining the renormalized 
coupling a = a(/i). We evaluate eq. (|3.33| ) in two schemes: The MS scheme, in which 



case Greg{u) is an entire function |T^, and the so-called MOM scheme, where the 



renormalized coupling is defined by 



(3.34) 



where R denotes any renormalization scheme or bare quantities. In this scheme we have 
F{u) = G{u) and eq. (|3.33| ) takes the particularly simple form 
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B 



MOM 



PI 



dv V 



u — V 



(3.35) 



Contrary to the MS scheme, the finite subtractions are themselves factorially divergent 
in this scheme. 

The right hand side of the first fine of eq. (|3.33| ) is easily evaluated as m — 1, given 
that 5rjS](f ) ~ 1/(1 - v) close to v = 1. We get 



§ /;7^(^rjS](^)-i?rjS](.))"=^-^i?rjS](n)(ln(l-n) + l) (3.36) 



Pi 



up to C((l — M)ln(l — u)) terms. To evaluate the second line in eq. ( p.33|) , we note 
first that we can drop the contribution from G^eg as long as we are not interested in 
singularities weaker than 1/(1 — u). In the MS scheme, this follows from Gj.cg{u) being 
entire, so that one gets only a ln(l — -u) -singularity from v close to one. In the MOM 
scheme, we also observe singular behaviour when v is close to zero, but since Bf^[T,]{v) 
is finite as v — >■ 0, one gets again only ln(l — u). Therefore, to the present accuracy, the 
results in the MS and MOM scheme are the same. Next we separate the double pole of 
Freg{u) at n = 1, since it gives rise to a logarithmic singularity of the f -integral at u = 1. 



In the remaining terms we can set u 
we obtain 



1. Adding the result from the previous equation 



5. [S](u)-=^i?rjS](n) 



Y Pi 




A PI 





"■prop 



(3.37) 



where 



6K, 



prop 



2 233 , 
-- + In 2 

3 12 



157 
I2" 



ln3 



-24E 

A:=3 

-1.7266 



(-l)^fc 

{k^ - 1)2 



1 



1 P 

2/fc2 + 2 fc2 



1 3A;2 
-- + 



arctan- 



1 



arctan- 



(3.38) 



Recall that /3i//5^ = 9/{ANf) = 0{1/Nf). We can now combine eqs. (^^), ( p30D and 
(|3.37|) to obtain the final result for the self-energy to subleading order in l/Nf. 



B[n{u) 



47r 



^2 



1 1 



21-u 



c+d ^ "-"-prop 



.oUi-u)Hi-u),j^ 



(3.39) 
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Figure 4: Cancellation of leading singularities between vertex and propagator correc- 
tions. Since these graphs should be inserted as subgraphs, the external photon lines are 
drawn as chains. 



with SK^l^ given in eq. ( p.31|) . By comparison of eq. ( p.ll| ) with eq. ( p. 101) , this result 



can be translated into the large-order behaviour of the self-energy. Notice that there is a 
cancellation of ln(l — u)/{l — u) terms between diagrams d and e of Fig. |^ and the only 
remaining logarithmic enhancement is proportional to /3i//3o- This remaining term can 
be traced to eq. ( p.36|) and then further to the coefficient of the pole of F{u) at u = 0. 
It is therefore unambiguously identified as the expected correction to b in eq. ( |1.2D due 
to two-loop evolution of the coupling (see Sect. 5). 

The cancellation between vertex and propagator corrections is of general nature and 
recurs in further examples in Sect. 4. It can be phrased as a cancellation between the 
two subgraphs shown in Fig. ^. Note that these subgraphs are of the same order in 
and for any diagram that contains the left subgraph there is a diagram that contains 
the right one and vice versa. Assume that these diagrams have N and -|- 1 chains 
respectively. The contraction of the left vertex subgraph gives c/(l — u^) with some 
constant c. Now consider the other subgraph. Contraction of the vertex subgraph in the 
inner box gives c/(l — un+i) with the same constant c when un+i is close to one. The 
subsequent contraction of the second box turns the simple pole into a double pole. If one 
takes into account all factors, including the different power of Po in eq. (|2.13| ), one finds 



that the coefficient of the double pole is {—i)c. After this contraction, the lower chain 
carries regularization parameter un + un+i- A factor (— z) from one photon propagator 
is obtained from joining the two chains into one (this can be done, since both have the 
same momentum). One can then make use of the (5-function constraint in eq. ( p.l3| ). 
From the left diagram we get simply 

(3.40) 



1 - (n - . . . - UN-i) 
and from the right 



It— til — ItjV — 1 







duNjz '-^ ^. (3.41) 

(1 - (U - Ul . . . - UN-l - UN)r 



Both contributions cancel each other up to non-singular terms. As a result the leading 
singularity drops out in the sum of the two diagrams that contain the graphs of Fig. 11 
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as subgraphs. This cancellation does not apply to the next-to-leading singularity of the 
diagrams. 

Finally we note that eq. ( |3.33| ) is universal and can be used for any Green function 
with the obvious replacement for the corresponding leading order Borel transform. 

4 Vertex function and vacuum polarization 

In this section we turn to the classical example of the photon vacuum polarization. For 
reasons that will become clear it is useful to consider first the fermion-photon vertex 
function. We restrict ourselves to the leading singularity at u = 1. 

The photon vacuum polarization is gauge- independent. One can introduce a gauge 
parameter ^, so that the Lorentz structure in eq. (|2.15|) takes the familiar form for 
covariant gauges. The result must be independent of ^. In particular, one can take 
Feynman gauge and considerably simplify the algebra connected with numerators. The 
price to pay for this simplification is the presence of singular factors like 1/((1 — u)uj) in 
individual diagrams, which are absent in Landau gauge, because the one-loop self-energy 
and vertex function are separately UV finite in this gauge. To render each diagram 
separately finite, one would have to add the usual counterterms. In practice we do not 
have to bother about these subtractions, since we know that all such singular factors must 
cancel in the final result. We have done the calculation for the vacuum polarization in 
an arbitrary covariant gauge and verified independence of ^ explicitly. Gauge-dependent 
intermediate expressions that follow below are given in Landau gauge, as before. 

4.1 Vertex function 

The result for the leading order diagram with one chain. Fig. ^b, is obtained from eq. ( p.5| ) 
after restoration of overall factors: 



BLo[r]iu 



l-u 



+ \ (if +p''^) liJi+iivii^ 



(4.1) 



The corresponding large-order behaviour is proportional to /^qTz!. Recall that the coup- 
ling i^—ig) to the external photon is not included in the Borel transform. 



4.1.1 Vertex diagrams with two chains 

The diagrams with two chains are shown in Fig. |[ Any maximal forest of any of these 
diagrams has two IPI elements, one of which is the diagram itself. According to eq. ( |3.2| ), 
the leading singularity arises from the combinations 



1 1 



(1 — ti)^ ' 1 — n 1 — ' 1 — uuj 



(4.2) 



where the last combination is in fact absent. The first combination arises from 'box' 
subgraphs as discussed in Sect. 2.2. For each diagram, in each forest, there are three 
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Figure 5: Vertex diagrams with two chains. 



choices for the second IPI element 7 7^ F, which is a one-loop subgraph. For each of 
the 18 distinct combinations we have to consider the region 14(7) = 1, where ^(7) is the 
sum of M-parameters of the lines of the corresponding one-loop subgraph. The result is 
denoted by (Nx), where N=l, . . . , 6 for the six diagrams and x=a,b,c for the three pos- 
sible one-loop subgraphs. Seven of the 18 combinations lead to insertions into reduced 
graphs that vanish. The remaining eleven ones are described as follows, according to 
their subgraph 7: 

(la): The 'small' vertex consisting of lines with momenta k2,p' — ki — k2,p — ki — k2- 
(lb): The 'box' subgraph. 

(2a) : The subgraph consisting of lines with momenta k2,p — ki — k2,p' — ki — k2,p' — k2. 
(2b): The subgraph symmetric to (2a) with lines k2,p — ki,p — ki — k2,p' — ki — k2- 
(2c): The 'crossed box' subgraph. 

(3a): The subgraph with lines ki,p — ki — k2,p — ki,p' — ki. 
(3b): The vertex subgraph on the upper fermion leg. 

(4a): The subgraph with lines ki,p — ki,p' — ki,p' — ki — k2, symmetric to (3a). 
(4b): The vertex subgraph on the lower fermion leg. 
(5a): The self-energy subgraph on the upper leg. 
(6a): The self-energy subgraph on the lower leg. 

In addition to subgraphs of self-energy [(5a), (6a)] and vertex type [(la), (3b), (4b)], which 
have superficial degree of divergence a; (7) = 0, we also encounter subgraphs with four 
external fermion legs [(lb), (2c)] and two fermion and two photon legs [(2a), (2b), (3a), (4a)] 
with c<j(7) = —2. 

Specifically, for the box-type subgraphs one obtains 



3tt 



1-uNf 



3 3 
±4 7m ® 7'' - 2>T5 ^ ^^^^ 



(4.3) 
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with the upper sign for (lb) and the lower sign for (2c). In the sum of both contribu- 
tions only the structure 7^75 ® 7'^75 survives. Similar cancellations occur between other 
subgraphs and one should combine them before insertion in the reduced diagram. It is 
then straightforward to calculate the combinations 



(I) 


= (la) 










(11) 


= [(2a)+(3a)] + [(2b)+(4a)] 






1 11 


/Pm) - 2p • p'7^ - ^ ^7^ / 


(III) 


= (3b)+(4b) = ^(QV-0QA<) 


(4.4) 




-\ip^ + p'^) > + ^ (^?'m+ /Pm) - \ 




(IV) 


= (5a)+(6a) 






1 13 

= Y2 ^'^^ ~ 3 ^^^'"^ 2 ^ 


/p^) - 2p • p'7^- ;^7^ / 


(V) 


= (lb)+(2c) = 2(g2^^-0Q^) 





It is understood that each line is multiplied by l/(167r^) and 1/(1 — u)(\ — u\) except 
for (V), where the latter factor is replaced by 1/(1 — uf' . Note that (I) and (IV) cancel 
against each other. The other terms add to 



(l)-(6) 



[r](n) 



LS 



1 



1 



1 



127r2 \\ — u \ — u\ 
3 1 



2 (l-n)2 



(4.5) 



The first line has exactly the same structure as the leading order result, eq. ( [4.1|) . After 



integration over the -u-parameters, it leads to a singularity of the form ln(l — ■u)/(l — -u), 

B^omn)^\^{\-n). (4.6) 

This is not the only contribution of this form, however. As in Sect. 3.5, eq. (|3.37|) , the 
universal propagator corrections add 



«) , (4.7) 



and by the same cancellation as for the self-energy only the term with I3i/ (3q is left. 

Adding everything together, we find that, including diagrams with two chains and 
propagator corrections, eq. ( ^?T|) is extended to 
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i?NLo[r](n) -j^2^''{Y^{^-^^Hi-n))\-l{Qh,-iQQ,) + l{{p'+p 



^ ^ (.y 'V,. — u,/i^/ ) -i- — \ \ 1) -+- r' 



7m 



The first line has a similar structure as in case of the self-energy. The two-loop evolu- 
tion of the coupling leads to a logarithmic enhancement in the large-order behaviour, 
PquI Inn, while being formally of order 1/Nf. The second line incorporates the con- 
tributions from box type subgraphs, which were absent for the self-energy (at the level 
of two chains). While this contribution starts at order 1/Nf, it leads to an enhanced 
divergence in large orders, PquI n for the reasons explained in Sect. 2.2. In the language 
of operator insertions, introduced in Sect. 1 and pursued in Sect. 5, we can identify 
this contribution with an insertion of the operator V^7/x75V''^7'^75'^ place of the box 
subgraph, see eq. (|4.^). The pattern of large-order behaviour for the vertex function 



with two chains is the same as for the vacuum polarization in |jTj] . With the help of IR 
rearrangement one can also compute the NLS, i.e. all terms that yield — u). In 
particular, the contribution from propagator corrections SKp^-op is the same as for the 
self-energy in Sect. 3.5. 

4.1.2 Vertex diagrams with one fermion loop and three chains 

In this subsection we analyze the set of graphs, labelled (1) to (6) in Fig. These 
diagrams have three chains, but contribute at order because the fermion loop 

brings a factor Nf. This set of diagrams is separately gauge- independent and we use 
Feynman gauge to simplify the algebra. The leading singularity is determined by the 
contributions with the maximal number of singular factors in eq. ( |2.23| ) which is three for 



the diagrams in Fig. ^ rather than two as for those in Fig. |. Because of fermion loops, 
the counting in and chains is no longer equivalent beyond a certain order. In general 
it is the number of chains and not the power in 1/Nf, which determines the strength of 
the renormalon singularity and the divergence in large orders of perturbation theory, so 
that the diagrams in Fig. ^ can be expected to dominate the large-order behaviour at 
order We determine the singularity at u = 1 up to terms ln(l — n)/(l — u). This 

includes those parts of the next-to-leading singularity that produce 1/(1 — u)^. 

As usual, the leading singularity is determined by subsequent contractions of one-loop 
subgraphs. The eight subgraphs that contribute to the singularity at m = 1 are enumer- 
ated (la) to (6b) in Fig. ^. There are many more subgraphs, which after contraction 
lead to vanishing (tadpole-type) reduced graphs. Other subgraphs like those in the last 
row of Fig. ^ have non- vanishing reduced graphs, but do not contribute at m = 1. For 
the first two graphs in that row, the contraction can be seen to correspond to operators 
of dimension eight, which contribute to m = 2 and higher poles only. In other words, the 
degree of divergence is uji^j) < —3. The last graph need not be considered, because the 
subgraph contains only fermion lines and no -u-parameter. 

The subgraphs of all eight non-vanishing contributions have the form of boxes or 
'crossed boxes'. When each box is combined with its corresponding crossed box, con- 
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(1) 

(5) 



(la) 



(4a) 




< 




< 




(2) 



(6) 



(3) 



(Ife) 



(2a) 



(5a) 



(6a) 



<C[ Kit -41 




(4) 



-<s: <£! Kff <h: 



(3a) 



-41 -<m: <m: 



(66) 



Figure 6: (l)-(6): Diagrams with three chains that contribute to the vertex at order 
1/A^|. (la)-(6b): The eight subgraphs which contribute to the leading singularity at 
u = 1. The last line shows subgraphs with non- vanishing reduced graphs that do not 
contribute at m = 1. 



traction of these subgraphs reduces to the insertion of '4''Jfi'j5ipip'j'^'j5'4'^ see eq. (|4.3| ). The 
four reduced graphs are shown in Fig. 0, where the black box is accompanied by a factor 
1/(1 — Mj — Uj) refer to the labels of the two contracted chains). Each of the graphs 
(/) to {IV) has two subgraphs with non-vanishing reduced graphs, called 'type A' and 
'type B' in Fig. |^. The contraction of type A yields another factor 1/(1 — Ui — Uj) 
being the same as in the previous contraction), while the contraction of type B produces 
1/(1 — u). The final contraction of the one-loop graphs to the right of the arrows in 
Fig. 1^ gives 1/(1 — m) in both cases. After integrating over the Ui with 



duidu2du3 5{u — ui — U2 — U3) 



(4.9) 



the leading singularity at m = 1 is of form 
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{I) 



type A: 



{11) 



{III) 



(IV) 




type B: 





Figure 7: (/)-(/ V): The reduced graphs corresponding to contraction of the subgraphs 
in (la)-(6b) in the previous figure. Each reduced graph has a subgraph of type A and 
B, whose contraction gives rise to the reduced graph shown on the right side. 



type A : 
type B : 



(l-tx)2 

ln(l - u) 

(1-U)2 



(4.10) 



In fact, we find that for each of the four diagrams of type A the traces vanish and there 
is no double pole contribution from this set. The contraction of type B results only in 
the tensor structure 7^ ® 7^^ but no 7^75 7'^75. As evident from Fig. ^ the pole term 
obtained in this contraction is related to operator mixing between ip'j^'y^ipip'j^'j^ip and 
tp'jfj^'ilj'ilj'-f'^ip. The final result for the leading singularity is then 



27 In(l-ti) 



47r^2 2Nf (1 



{Q\- iQQ^) , 



(4.11) 



where Q is the momentum of the external photon. 

To obtain all terms 1/(1 — m)^ in the final result, we have to collect in addition all 
contributions with two singular denominator factors 1/(1 — u), since the contributions 
from terms with three singular factors (type A above) vanish. Since one such factor 
comes from the last contraction, all such terms arise from two-loop four-fermion graphs, 
obtained by cutting the fermion loop in the diagrams (1) to (6) in Fig. The problem 
therefore reduces to computing first the leading and next-to-leading singularity for the 
four-fermion Green function with three chains, extending eq. (|4.3| ) to the next order in 
the chain expansion. The resulting expression is then inserted into a reduced graph of 
type B, shown to the right of the arrow in Fig. |^. 
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As in the case of the self-energy, we use the method of infrared rearrangement to 
obtain the coefficient of the singular factor 1/(1 — m) of the two- loop (three chain) four- 
fermion graphs. Since the singularity is local and proportional to Q^, no differentiation is 
required. Therefore we can set all external momenta to zero from the start and introduce 
a new external momentum p that flows through a single line. We then encounter integrals 
of the form 



(27r)4 (27r)4 {{p + A:i)2)3-«i (A;2)3-«2 ((/t^ + A^a)^)!- 



"3 



which are easily evaluated. Then, as a straightforward generalization of the formulae of 
Sect. 3.4, one subtracts the contributions from the leading singularity before integration 
over M-parameters. We can set m = Ui + ^2 + = 1 after this step and perform the 
integration. We then obtain for the Borel transform of the two-loop four-fermion graphs 
near u = 1 

i Stt 97 

[ln(l-n)+Cbox]7M®7^ (4.13) 



l-uNf 2Nf 



where Cbox = 1.016 . . . results from numerical integration over a product of F-functions. 

The corresponding contribution to the vertex follows from inserting the previous line 
in place of the black box in the last diagram of Fig. |^. The result is 



Br,,J^]{u) ^ [ln(l - u) + (Q%,- ^Q,) , (4.14) 



extending eq. ( [4.11|) including all contributions of order 1/(1 — u)"^. Note that these 



diagrams produce the dominant pole in the vertex function at order '^/Nj. 
4.1.3 Summary 



Collecting all contributions at order ^/Nj, eq. ( |4.8| ) and eq. (|4.14| ), the final result for 
the large-order behaviour of the vertex function reads 



00 



r(p,p» = (-iff)^r„a"+i (4.15) 

n=0 



with 



n— KX) 



^ e^/?o" nlj (1 + ^(Inn - V(l))) [ - ^(QS^" + ^ ((/ + p" 



7^. 



+ F7m; 



+ J^n (^-lnn + V'(2)+C7box + o(^n ) {Q^-^,)}- (4.16) 



The missing 0{l/nlnn) terms come from the uncalculated ln(l — m)/(1 — u) terms of 
the diagrams with three chains and one fermion loop. 
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4.2 Vacuum polarization 



With the resuhs for the vertex function, the resuhs for the vacuum polarization are 
immediate. A specialty of the vacuum polarization is that the skeleton diagrams corre- 
sponding to diagrams with one chain have two loops. This gives an additional factor n 
in large-orders as compared to the vertex already in lowest order. 



4.2.1 One chain 

The diagrams with a single chain are part of Fig. |1|. The diagram where the chain forms 
a self-energy subgraph does not contribute to the leading singularity. The reason is that 
the only forest that does not lead to tadpole reduced graphs is the one with the self- 
energy subgraph. If p is the external momentum of this subgraph, the self-energy close 
to M = 1 is proportional to ^, see eq. (P7^), and the insertion of this polynomial into 
the reduced graph eliminates any l/p^ from loop lines of the reduced graph. Thus the 
reduced graph is effectively a tadpole and vanishes. We are left with the first graph in 
the second line of Fig. 0. The contraction of the vertex subgraph leads to the insertion 
of eq. (f4.1|) , so that 



i?Lo[n](n) 




(4.17) 



Notice that the terms involving p and p' in eq. ( [4.1|) do not contribute, when inserted 
into the reduced diagram, because they eliminate one of the two denominators. The 
overall factor two accounts for the two possible insertions into the left or right vertex. 
The second hne above reproduces| the asymptotic behaviour of the exact result of p. 
The corresponding large-n behaviour is proportional to jS^nln. To compute the 1/n 
corrections one would use IR rearrangement with additional subtractions of non-local 
terms in external momenta. 



4.2.2 Two chains 

We calculate the leading singularity for the diagrams with two chains, see Fig. ^ A 
maximal forest contains three IPI elements. Now we encounter for the first time forests, 
where non-trivial elements are disjoint and do not form a nested sequence, for instance 
for the third diagram in class (IV). However, close to m = 1, such forests do not contribute 
to the leading singularity, since they produce at most ln(l — u)/{l — u). When the non- 
trivial elements form a nested sequence, the reduced diagram is a tadpole, whenever the 
smallest subgraph in this sequence contains one or both of the external vertices. 

^ There is a factor ~Nf/{A'K) difference due to different conventions. Here the two coupUngs {—ig) 
to the external photon are not included by definition, while in Q the Borel transform of Il/a is given 
with a = aNf. 
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Figure 8: Vacuum polarization diagrams with two chains. The crosses denote insertions 
of external photon legs. 

With these restrictions in mind, let us turn to the diagrams (I) in Fig. |. The 
smallest graph in each non-vanishing maximal forest is a self-energy subgraph. Close 
to Ui = 1, where Ui is the regularization parameter of this subgraph, it reduces to an 
insertion proportional to ^, with p the external momentum of the subgraph. Again 
this insertion eliminates one denominator, so that the reduced graph in fact vanishes. 
One can convince oneself that for class (II) a non-tadpole reduced graph always requires 
contraction of the two-loop self-energy subgraph. We then get the same cancellation 
as for (I). Thus there is no contribution to the leading singularity from the first seven 
diagrams. 

Considering class (III) and (IV), we can divide all maximal forests into three groups: 
(i) A forest of type {^Vertex, r}, where J^vertex is a forest that has already been evaluated 
for the two-chain vertex (Fig. ^ and F is the entire diagram, (ii) Forests with nested 
IPI elements, not belonging to group (i). (iii) Forests with disjoint one- loop elements. 
It turns out that (ii) always leads to reduced graphs of tadpole type and (iii) does 
not contribute to the leading singularity as noted above. Thus we arrive at the same 
conclusion as in the case of a single chain that the leading singularity of the vacuum 
polarization is given by insertion of the result for the vertex function, in the present 
situation eq. ( |4.5|) or eq. ( |4.8| ), if one includes already the propagator corrections. Again, 



only the terms involving Q contribute for the same reason as above. Since such terms 
do not arise from class (III) diagrams, only the final four diagrams in Fig. ^ contribute 
to the leading singularity. (This could be deduced without knowing the result for the 
vertex function, applying the same argument as to class (I) and (II).) The insertion leads 
to exactly the same integral for the reduced graph as in eq. (|4.17|) . Therefore we get 
(including propagator corrections) 



l-§ln(l-u) 
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(4.18) 



and the corresponding large-order behaviour of the perturbative expansion (including 
now the 'light-by-light' scattering contributions) 



N 
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|l(lnn-V(2)) 
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-n 



8Nf 8Nf 



n{-lnn + V'(3) + C^ox} 



(4.19) 



The corrections to this result are 0{1/Nf ,111^ n/Nj) relative to the unity in brackets. 
The pattern of divergence is the same as for the vertex function in that contractions 
that lead to four-fermion operator insertions are enhanced by a factor of n as discussed 
in Sect. 2.2. In analogy with the vertex function, Sect. 4.1.2, the 'light-by- light' diagrams 
with three chains and two fermion loops, not shown in Fig. ^, contribute to the vacuum 
polarization at subleading order in 1/iV/. 

Setting the expression in curly brackets to zero, eq. ( |4.19| ) agrees with the correspond- 
ing result in the erratum to |jT^, where the light-by- light contributions are eliminated 



by a certain charge assignment to the fermions. Compared to the approach taken there, 
based on the background field technique, the diagrammatic formalism presented here 
lacks the elegance of maintaining explicit gauge invariance. This leads to a proliferation 
of terms such as in eq. ( f4.4| ), most of which cancel in the final result. On the other hand, 
the conceptual framework developed in Sect. 3 allows us to systematically compute cor- 
rections, either from more then two chains or contributions subleading in n for large n, 
such as the first correction to the normalization K of the large-order behaviour (see the 
example of the self-energy in Sect. 3) and thus elucidates the origins and organization of 
divergence for arbitrarily complicated diagrams. The understanding of the systematics 
of the chain expansion to all orders is the subject of the following section. 



5 UV renormalon counterterms and 
renormalization group 

The previous examples illustrate that diagrams with a larger number of chains (sup- 
pressed by powers of 1/Nf) lead to stronger divergence in large orders of perturbation 
theory. Except for an enhancement of n due to the new class of box-type subgraphs, 
the enhancement is logarithmic in n. For each additional chain one can obtain at most 
one Inn, as there is one additional singular factor in the w-parameters. The expansion 
parameter of the chain expansion is in fact \nn/Nf and any finite order approximation 
in 1/A^/ does not provide the correct asymptotic behaviour in n. 

The situation is reminiscent of standard applications of the renormalization group, 
which allow us to deduce the large-momentum behaviour of Green functions, which 
can not be determined from finite-order perturbative expansions in a. The key idea 
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is to prove factorization of the dependence on the variable under consideration and 
relate this dependence to one on an artificial factorization scale. In the present case, 
factorial divergence in n is equivalent to poles in the Borel parameter m, which in turn 
are related to the regions of large loop momenta. Specifically, for the contribution from 
a given forest, the leading singularity could be determined from the d^fc/fc^-piece in 
the expansion of the smallest non-trivial element of the forest in its external momenta 
and thereby associated with insertion of a dimension six operator. Generalizing this 
observation, the desired factorization would be summarized by the statement |Q that 
the leading UV renormalon could be compensated in all Green functions by adding 

4? = -iE^^(«(^))o- (5-1) 

to the renormalized Lagrangian, where Oi are dimension six operators and Ei{a{fi)) 
factorially divergent series in a with finite coefficients that depend on the renormalization 
conditions implied for the (usual) renormalized Lagrangian as well as the operators Oi. 

Such a statement implies that UV renormalons are local. Indeed, all explicit results 
for large-n behaviour given so far were straightforwardly obtained as local quantities 
(polynomial in external momentum). But the expectation to get each time a local 
result is too naive and logarithms of the external momenta arise at the level of 1/n- 
corrections.0 The situation is rather similar to the calculation of UV counterterms. If one 
does not remove sub divergences and calculates the pole part of a diagram in e = (4 — (i)/2 
(where d is the space-time dimension) the result is inevitably non-local. To get a proper 
answer one should first take care of all the subdivergences of the diagram by inserting 
the corresponding counterterms which subtract the subdivergences. In the case of UV 
renormalons, one should also subtract, for subgraphs of a given graph, the factorial 
divergence associated with subgraphs. This procedure is of course implicit in writing 



eq. ( |5.1| ) and the non-trivial assertion is that the remaining 'overall' UV renormalon 
divergence for the given graph is indeed local. The UV renormalon subtractions will 
be organized in the same way as usual subtractions of UV divergences that enter the 
standard UV i?-operation (i.e. renormalization at the diagrammatical level). 

In Sect. 5.1 we illustrate the procedure on the simplest example, then defining the 
general 7?.-operation in Sect. 5.2. The analogy with the usual i?-operation suggests 



eq. (5J-), but we do not actually prove this result. The resulting combinatorial struc- 
ture of the 7?.-operation enables us to straightforwardly translate our result into the 
Lagrangian language. The best way to establish this property manifestly is to apply 



formulae of the so-called counterterm technique PUI BD], which are essentially based 



on combinatorial properties of the i?-operation and locality of counterterms. We shall 
list the corresponding formulae in Sect. 5.3. These formulae do not explicitly give the 
operators which are added to the Lagrangian in every concrete situation. In Sect. 5.4 
we characterize these operators, the corresponding renormalon counterterms and derive 
renormalization group equations for the coefficient functions Ei that enter eq. ( |5.1| ). The 
solutions to these equations sum all contributions of order [Inn/Nf)'' and provide the 
correct asymptotic behaviour in n, up to an overall constant, that can be determined 



'See the discussion in Sect. 3.3 after eq. (3.13). 
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Figure 9: Leading order vacuum polarization: Diagrams (a) and counterterm insertions 
(b), (c). 

only order by order in 1/Nf. For the vacuum polarization, this solution is found in Sect. 
5.5. 

5.1 Counterterms: An example 



The simplest case of subtraction of sub divergences, considered already in |12], is the 
photon vacuum polarization with a single chain. 

Let Gj^{q] a) represent the perturbation series in a, generated by the set of diagrams 
r, such as shown in Fig. ^a for the vacuum polarization. Assuming that we know how 
to extract the factorial divergence for a given set of diagrams including 1 / ra^-corrections 
to the leading asymptotic behaviour up to some specified value /cq of ^! we define an 
operation K, that does this. If the factorial divergence comes from a branch point in the 
Borel plane as is generally the case, ko cuts off the infinite series of corrections to the 
leading asymptotic behaviour. If it comes from a pole, as for a single chain, this series 
terminates by itself. Since there is only one regularization parameter u, the singularities 
at M — *■ 1 are of the form 1/(1 — u)^ and 1/(1 — u). The residue of the simple pole 
of the corresponding skeleton diagram with analytically regularized photon propagator 
involves a logarithm of the external momentum. In terms of large-ra behaviour, we find 



X <'n-ln I -^e^' 1 + , (5.2) 




see eq. ( |4.17] ). The 1/n-corrections are taken from the exact result |p. Because of the 

ation can not be obtained from the insertion of 

^E'^{a{ii)) , 0'^ = d^F-Pd^F^p (5.3) 



logarithm, the previous equation can not be obtained from the insertion of 

1 



into the two-point function.^ 

To arrive at a local result for the overall counterterm associated with the vacuum 
polarization diagrams one should insert the (renormalon) counterterms A(7j) for two 

^^F^i, is the field strength tensor of the photon. We use 'primes' on coefficient functions and operators 
here, because later we introduce another basis of operators. 
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1-loop vertex subgraphs 71 and 72 of the first diagram in Fig. ^ and deal with the 
quantity 



7^'^^(<7; u) = {l + A(7i) + A(72))nr(g; n) 



(5.4) 



Here TZ' is the incomplete ('renormalon') 7?.-operation, i.e. without the last (overall) 
counterterm. The various symbols that appear in the previous equation will be defined 
more precisely in the following subsection. 

The counterterms connected with the vertex subgraphs can be deduced from eq. ( [4.1|) . 
Adding the term 



QlTfJ.' 



(5.5) 



to the Lagrangian, the vertex function F is free from the first UV renormalon to the 
leading order in 1/A^/.|^ This new term in the Lagrangian contributes to the vacuum 
polarization at leading order through the diagrams shown in Fig. Pd, which correspond 
to the two additional terms in eq. (^.4]). Their evaluation gives 



(A(7i) + A(72))nr((z;a) 



-V. 



AT, 



7_E + In 47r 



(5.6) 



The insertion of counterterms for subdiagrams generates new UV divergences. Their 
subtraction corresponds to a renormalization prescription for the set of dimension six 
operators. Since the w-parameters are not sufficient to regularize these UV divergences, 
we have applied dimensional regularization. Choosing the MS subtraction scheme, op- 
erator renormalization on the diagrammatic level is accomplished, when the insertion 
of the (finite) renormalon counterterm is followed by the action of the MS -R-operation. 
Thus we apply the product 



Rn' = l + {1 + A(r/7i))A(7i) + (1 + A(r/72))A(72 



(5.7) 



instead of eq. (^.4]). Here 1 + A(F/7i) = i?(F/7j), i = 1,2 are /^-operations for the 
reduced diagrams. Combining eqs. (|5.2|) and ( ^.6|) , we obtain 



/Ci^7^'^^(Q;a) 



1^' 



n + --C 
b 



Nc 



(5.^ 



"'^"'^To be precise, additional operator structures must be added to compensate the terms involving p 
and p' in eq. (4J). We may ignore these, because their insertion leads to vanishing reduced graphs in 



the present example, see also the discussion in Sect. 4.2.2. 
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(a) 



(b) 



Figure 10: A self-energy diagram and its counterterms. 



The logarithm of has dropped out and the overall counterterm, represented by dia- 
gram c of Fig. 1^, is local as required. It is straightforward to read off the coefficient E'^ 
from the previous equation, so that the vacuum polarization is free from the first UV 
renormalon, if the operators and O'^ are added to the Lagrangian with appropriately 
chosen coefficients. 

For diagrams with two chains, the diagrams with counterterm insertions also have 
to be calculated in the chain expansion. Consider the diagram shown in Fig. p!0| a that 
contributes to the self-energy. Combining results from Sect. 3, the Borel transform of 
this (set of) diagram(s) in the vicinity of m = 1 can be represented as 




(5.9) 



where SK^'^^ has been given in Sect. 3.4 and the unspecified constant is scheme- dependent. 
The two diagrams with insertions of O'^, Fig. contain one chain and its contribution 
is given by the product of the series for E2{a{fi)) with the series for the renormalized 
diagram with operator insertion. The Borel transform is given by the convolution 



2- — / duidu2S{u — ui 
Po J 





47r 



02 \ 

-^e^j P{u2)-R{u2) 



, (5.10) 



where both, P{u2) and R{u2), behave as —1/(2^2) as 1*2 — >■ and the factor of two 
accounts for the two identical diagrams. The function P{u2) arises from straightforward 
calculation of these diagrams with analytically regularized photon propagator. Since the 
diagrams are quadratically divergent, the series of ultraviolet poles of P{u2) starts from 
U2 = —1,0,1,.... The function R{u2) incorporates the ultraviolet subtractions in the 
adopted scheme for the operator renormalization (in this case mixing of into ipd'^^i/j). 
In the MS scheme only the logarithmic divergences are subtracted, so that R{u2) is 
analytic except at U2 = 0. The function R{u2) can be calculated with the methods of 
jlT|. Using the series for E'.^ given above, the previous integral can be evaluated and 
yields 




1 1 
2'Wt 



ln(l 



In 



>2 



+ const' 



(5.11) 
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The sum of eqs. ( ^.91 ) and can again be absorbed into a local counterterm pro- 



portional to '4)8'^ ^ip. The logarithm of arises only from Ui 1, and the unspecified 
constant depends on the subtraction scheme for the dimension six operators. 

It appears that there is also a logarithmic singularity in u from the other integration 
boundary, U2 1, because P{u2) is singular at this point. This contribution would be 
proportional to QiQ"^ I IJ^Y ■ However, we can always choose the series for £'3 to start at 
some sufficiently large no, so that the region M2 — >■ 1 is suppressed by a zero of B[E'.^{ui) 
and is non-singular. 

It is worth emphasizing that the non-local terms are related to logarithmic divergences 
of the reduced diagrams (diagrams with insertion of dimension six operators), while 
their coefficient is also unambiguously related to the leading singularity of the original 
diagram. (As in dimensional regularization a pole in e is accompanied by a logarithm 
in momenta.) This property allows us to sum the leading singularities to all orders in 
with a variant of standard renormalization group techniques. 

We also note that as long as we consider only the first UV renormalon at u = 1, we 
do not have to consider multiple insertions of dimension six operators. For instance, for 
two insertions, to be sensitive at once, to the singularity at u(7i) = 1 and at ^(72) = 1, 
we need at least m = 2, since the 71 and 72 can not have common u-parameters. The 
same conclusion follows from 



2n 



(/?o"n!) (/3o"n!) ~ 2 (^^ j {2n)\^ (5.12) 
in the language of large-n behaviour. 
5.2 7?-- operation 

The illustrated counterterm procedure can be generalized to arbitrary diagrams. More- 
over, since according to eq. (|2.23| ) the combinatorial structure of poles in the w-parameters 
that eventually give rise to singularities in u of the Borel transform is identical to the 
combinatorial structure of logarithmic UV divergences, the 'renormalon 7?.-operation' 
can be constructed in analogy with the usual /^-operation. 

Let us recall that the standard i?-operation applied to a Feynman integral Fr acts 
31], |2|, H, g 



as 



RFr= A(7i)...A(7,)Fr = i?'Fr + A(r)Fr, (5.13) 

where A (7) is the corresponding counterterm operation, and the sum is over all sets 
{71,..., 7j} of disjoint divergent IPI subgraphs, with A(0) = 1. The 'incomplete' R- 
operation R' by definition includes all the counterterms except the overall one A(r). 
Within dimensional renormalization, the action of the counterterm operation on the 
given Feynman integral is described as 

A(7) Fr = Fr/^ o P^, (5.14) 
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where -Fr/7 is the Feynman integral corresponding to the reduced graph r/7 and the 
right-hand side of eq. ( ^.141) denotes the Feynman integral that differs from -Fr/7 by 
insertion of the polynomial P^yiq,]]!) into the vertex to which the subgraph 7 was 
collapsed. In the MS-scheme the coefficients of these polynomials are represented as 
linear combinations of poles in e = {4—d)/2. It is implied that dimensional regularization 
is still present in eqs. (|5.13|) and ( |5.14|) . 

In the framework of the renormalization schemes based on dimensional regularization 
the Pjigyin) are polynomials with respect to masses m and external momenta g of 7 
(remember that it is sufficient to consider the pure massless case in our problem). 
The degree of each equals the degree of divergence 0^(7). In the MS scheme these 
polynomials are defined recursively by equations 

= A(7) = -k^R' (5.15) 



for the graphs 7 of a given theory. Here K^^ is the operator that picks up the pole 
part of the Laurent series in e. Note that the essential part of the basic theorem on 



the /^-operation in the framework of dimensional renormalization pO[ is just the above 
polynomial dependence of diagrammatic counterterms P^yiq,]]!) on masses and external 
momenta. 

Let us now define an operation TZ, with the same combinatorial structure as R in 
eq. ( ^.13|) , which removes a (finite) number of terms in the large order behavior of a 



given set of diagrams (that is, the leading large-n behaviour, including 1/n'^-corrections 
up to some ko, is removed). This operation is applied to sets of diagrams T, defined as in 
Sect. 2.1 except that the photon propagators in T are chains rather than complete photon 
propagators. We assume that this set of diagrams includes the usual UV counterterms, if 
necessary, so that Gfiq] a) is a series in a with finite coefficients, when the regularization 
is removed. We define the TZ recursively by 

7^G^= Yl ~^{li)---~^{lj)Gf = n'Gf+~K{G)Gf, (5.16) 

7iv,7j 

Hl)Gf = ^G^oV^, (5.17) 

where the sum is over all sets {71, . . . , 7^} of disjoint IPI skeleton subgraphs with degree 
of divergence 07(74) > —2. Every polynomial is represented as a finite sum of terms 
of the form 

Y.a^Ponlnnn^npjiiq') (5.18) 

n 

where pji are polynomials of degree 0^(7) + 2 in external momenta of 7. These countert- 
erms P-y are determined by 

= A(7) = -KRn'G^{q, a) . (5.19) 
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The operator K extracts the factorial divergence (including l/n'^-corrections up to the 
chosen ko) of the series in a to its right. The incomplete 7?.-operation TZ' is described by 
eq. ( |5.16| ). It is implied that the initial /^-operation in the MS-scheme that subtracts the 
(usual) UV divergences in 7 is implicit in G^{q,a). (Alternatively it could be included 
in the definition of IZ. In terms of Borel transforms these subtractions correspond to 
the subtraction of poles at ^(7) = plus a series in 'u(7).) The insertion of V for 
subgraphs of a given graph generated new UV divergences, which are subtracted by 
subsequent application of i? = i?'^^, the usual renormalization in the MS-scheme. Note 
that these extra MS-counterterms (associated with the renormalization of dimension six 
operators) are not inserted in the usual way: We have counterterms accompanied by 
such counterterms but do not have these counterterms alone. 

When calculating the counterterms with the help of eq. ( |5.19|) we introduce di- 



mensional regularization which will be switched off only after the action of the operation 
R but before extracting the large-order behaviour with /C. As in the usual calculation of 
counterterms IR rearrangement can be used (as in Sect. 3.3) to simplify the calculation 
of V^, i.e. we differentiate a;(7) + 2 times in the external momenta, then put them to 
zero for all resulting diagrams and introduce, in the simplest way for the calculation, a 
new external momentum. If there is no way to avoid IR divergences when following this 
prescription, we can apply the so-called i?*-operation which removes these spurious 
divergences. In this case, we have (after differentiation) the product R*TV in eq. (|5.19|) 
instead of RTZ' . Note that it is natural to apply the dimensional R* because we then get 
zero not only for massless vacuum diagrams and tadpoles themselves but as well for the 
corresponding 'i?*-normalized' values. 

We state without detailed proof that the will indeed be local (polynomial), so that 
the corresponding subtractions can be implemented in terms of dimension six operators 
in the Lagrangian. A complete proof would follow the same strategy as for the usual 
/^-operation. 

5.3 Counterterms in the Lagrangian 

The operation TZ was defined at the diagrammatic level. However, the combinatorial 
structure of the /^-operation, eq. ( ^.13| ) together with the locality of the counterterms 



enable one to express the renormalization procedure by inserting counterterms into the 
Lagrangian. An explicit way to do this is to apply formulae^ of the counterterm tech- 
nique |29, 2^. For the generating functional of Green functions the basic formula 



looks like 

i?exp{i5(0)} =exp{iS,(0)}, (5.20) 



where (f) stands for the collection of all the fields of a given theory, 5* is the interaction 
part of the action and is the renormalized action (interaction plus counterterms) 

^^For brevity, we list the formulae for normally ordered Lagrangians. In the case without normal 
ordering which is really implied, the corresponding generalizations differ by inserting some additional 
operator - see details in ^ . 
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which is exphcitly expressed through the counterterm operation as 



Sr{(t)) = A{exp{iS{(p)} - 1). (5.21) 

By definition, the -R-operation and counterterm operations, as apphed to functionals, 
act on each diagram involved in their diagrammatic expansions. 

Now we observe that the same combinatorial arguments that resulted in eqs. ( [5.20| ) 
and ( ^.21| ) can be used to describe the 7^-operation in eq. ( |5.16|) in the Lagrangian 



language. We have 

7^exp{^5(</))} = exp{z5^((^)}, = A(exp{i5(0)} - 1), (5.22) 

where application of the counterterm operation A to the functional in the right-hand 
side of eq. (|5.22|) reduces to the action of the diagrammatic counterterm operation (given 
by eqs. ( |5.17| ) and ( [5.19D ) on the whole classes of diagrams with vacuum polarization 
insertions into the photon propagator. 

Eqs. (|5.22| ) provide an explicit realization of the fact that the renormalons can be 
compensated by inserting counterterms into the Lagrangian. The calculation of these 
counterterms is performed by use of eq. (|5.19| ) - see examples in the previous section. 



Remember that to absorb the large order behavior into the Lagrangian we restrict our- 
selves to the dimension six finite extra counterterms while the infinite extra counterterms 
serve to renormalize diagrams which are generated by insertion of these dimension six 
operators. 

5.4 Renormalization group equations 

Since we do not necessarily want to modify the QED Lagrangian by higher-dimensional 
operators, the main relevance of the statement that UV renormalons could be com- 
pensated in this way is that it provides information about the large-n behaviour of 
perturbative expansions in the unmodified theory beyond the 1/A^j-expansion. 

To develop the argument it is convenient to adopt a different language for the coun- 
terterms (and we stress that the physical content stays the same - in particular it is 
entirely perturbative, although some expressions may not appear so). Given the Borel 
transform B[G]{u), we can formally recover G{a) by 

oo 

G{a) = / dne-"/(^°°) B[G]{u) . (5.23) 







Let us assume that the perturbative coefficients of G{a) are made real by multiplication 
of appropriate factors of i so that B[G]{u) is real. Because of the singularities on the 
integration contour, starting at m = 1, we define the integral with a contour in the 
complex plane slightly above the real axis. G{a) acquires an imaginary part, whose 
dependence on a is in one-to-one correspondence with the nature of the singularity. If 
B[G]{u) has a simple pole at u = 1, 
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B[G]{u) ~ lmG{a) = ^ g-^/^'^"") . (5.24) 

1 — u po 

A double pole gives an additional factor l/(/?oa)- Instead of choosing E'^ as series in a 
to compensate factorial divergence, let us now choose E'^ to compensate the imaginary 
parts of the Borel integral generated by factorial divergence. In this language E'^ in 
eq. ( p. 51) is replaced by 

-^^^:e-V(/5o"). (5.25) 

/^^ PO DTT 

By construction the Green functions computed from the modified Lagrangian (i.e. in- 
cluding dimension six counterterms) satisfy 

ImG'^°^(a) =o(e-2/('^o°)) , (5.26) 



where the remaining imaginary part is related to singularities at m = 2 and higher. Thus, 
ImG(a) = -i^^i?,(a)Go,(a) + 0(e-2/(A,a)^ (5 27) 



for the Green functions in the unmodified theory. Here Go. (a) denotes the Green func- 
tion G with one zero-momentum insertion of the dimension six operator Oi. As discussed 
in Sect. 5.1 multiple insertions contribute only to singularities at m = 2 or higher. A 
simple equation like the previous one can not be expected to hold for these singularities. 

A Green function G{a), and therefore ImG(a), and Go. (a) satisfies a standard renor- 
malization group equation. Comparing the renormalization group equations for ImG(a) 
and Go- (a), we obtain the evolution equation for the coefficient functions. 



IJ'' 



d_ 

da 



Sij + Ijiia) 



Ej{a{n)) 



0. 



(5.28) 



where 

HM = {Z-^)k, , or = Zi,0^ (5.29) 



are the anomalous dimension matrix and renormalization constants for the set of di- 
mension six operators. (We assume dimensional renormalization so that operators with 
different dimension do not mix.) Since Ei{a{fi)) does not depend on explicitly, we 
arrive at 

Ej{a) = 0, (5.30) 
where j'^ denotes the transpose of 7. The solution is given by 
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E,{a{n)) = exp / da ' ^. 7. Ej (5.31) 



and a-independent initial conditions Ej. Consequently the a-dependence of ImG'(a) 
and therefore the n-dependence of the large-n behaviour of the perturbative expansion of 
G{a) is completely determined by renormalization group considerations. (To obtain the 
complete l/n'^-corrections to the leading large-n behaviour one has to compute Go. (a) 
to order in addition.) The initial conditions Ej provide the overall normalization, 
which has to be determined, for each operator, from a suitably chosen Green function. 
This calculation can be done systematically only through expansion in 1/Nf. A solution 
of form ( |5.31| ) has been obtained previously by Parisi [Q and employed in |^ in a heavy 
quark effective theory context. Similar ideas, although in a technically somewhat differ- 



ent set-up, were used by Vainshtein and Zakharov to find the asymptotic behaviour 
of a current two-point function in a QED-like model. 

Let us now discuss the set of operators O,. Since we consider operator insertions into 
Green functions, which can have off-shell external momenta, we keep operators which 
vanish by equations of motion as well as gauge-variant operators in the general case. 
However, not all operators are of importance for any given Green function G, since the 
contribution to ImG from an operator Oi with coefficient Ei (which is independent of 
G) can have additional factors of a from Go. (a), which lead to 1/n-suppression in the 
contribution to the large-order behaviour of G. 

Specifically, in the example to follow in Sect. 5.5, we consider the operators 



Ci = V57/.75V'V^7''75V' (5.32) 
O2 = V'7m'0V'7^V' , 

C4 = g^D,F^'P{D^'F^p + gi;-fpi;) 



Note that and O4 vanish by the equation of motion for the photon field, but should 
be kept, because for the vertex function and vacuum polarization we consider an off- 
shell external photon. We omit the operators that vanish by the equation of motion 
for the fermion field from the above list, because in the examples to follow we assume 
the external fermion lines to be on-shell for simplification. We have also omitted the 
gauge-variant operators. 



5.5 Asymptotic behaviour 

5.5.1 Solution to the renormalization group equation 



We now solve the renormalization group equations, eq. ( |5.28| ), which allows us to derive 



a 'renormalization group improved' form for the large-order behaviour of the vertex 
function and vacuum polarization which replaces the next-to- leading l/iVy-result given 
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in Sect. 4. This improved form sums all (In n/A^j) '^-contributions to all orders in 
and provides the correct asymptotic behaviour in n with corrections being of order 1/n 
or 1/Nf. In the case of the vertex function, we consider only the case when the external 
fermions are on-shell, so that the Dirac structures involving p and p' in eq. (|4.1| ) can 
be neglected. As a consequence, eq. ( p.32| ) exhausts the hst of all possible dimension 
six operators, since we have to consider zero-momentum insertions only and restrict 
ourselves to gauge-invariant quantities. 

In the operator basis chosen above, the anomalous dimension matrix is block trian- 
gular, since the operators that vanish by the equation of motion do not mix into those 
that don't. The result for the anomalous dimension matrix is 
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0-a 


3a 
27r 
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2Nf+l a 
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3 vr 
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2Nf+l 


NfO 








~3^ 
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Nf 
127r^ 






Nfa 
3tt 



\ 



(5.33) 



All entries are given to lowest order in a for each entry. Corrections of order a will 
generate only l/ri-corrections to the asymptotic behaviour in n and we do not consider 
such corrections. Eq. (|5.3CI| ) is easily solved with this anomalous dimension matrix. It 
is possible and advantageous to choose the solution such that a appears only in the 
combination a = f3oa. We first solve for Ei and E2 by diagonalizing the 2 x 2-matrix 
pertaining to operators Oi and C2. Then we integrate to obtain first E3 and then E4. 
Introducing 



^1,2 



2Nf + 1 
6 




7i 



.A 

7r/3o ' 



(5.34) 



we obtain 



Ei{a) = e^i/^'a-^i/^o (1 + 0{a)) | 



f2Nj + l 



-C/ 1 H Eo 

11 



-1/2 



+ 11 



H Eo 

11 



(5.35) 



E2{a) = e-^/^a-'^i/'^o (1 + 0(a)) 



E-\ H Eo 

11 



'2Nf + 1 




E^ -\ Eo 

11 



(5.36) 
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Ei(a) 
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(5.38) 



£'1 H Eo 

11 



These results can be translated into the original language, where the Ei{a) were facto- 
rially divergent series, by the rule 



vr 



(5.39) 



One can then determine the large-order behaviour of any given Green function directly 
from 



G{a) Go,{a)- 



(5.40) 



5.5.2 Matching 

It remains to determine the integration constants Ei by matching the solution to the 
renormalization group equation above with the results of explicit calculations in the 
chain expansion. We do this interpreting Ei{a) as factorially divergent series in a. Com- 
parison with explicit calculation provides a non-trivial consistency check as it must be 
reproducible with four n -independent constants E^. 

Four-fermion scattering. Since the insertion of O3 and into the four-point 
function vanishes, the four-point function allows us to determine Ei and E2. As men- 
tioned before, we consider on-shell external fermion legs. Most of the ingredients for 
the calculation of four-fermion scattering to next-to-next-to-leading (NNLO) order in 
are scattered among previous results and we quote results only after putting them 
together. 

At leading order, we have a tree diagram with a single chain. It does not yield 
divergences in large orders. At next-to-leading order, the set of IPI diagrams consists of 
the box diagrams. From eq. ( [4. 3] ) we obtain 

^lPI,NLO n^oo _J_,C2 (^_^^ /3o"n!7M75 7^5, (5.41) 
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where r„ is the coefficient of za""*"^ in the perturbative expansion and the factor of two 
comes from an equal contribution of scattering and annihilation type diagrams. There 
are two contributions from one-particle reducible (IPR) diagrams. One from a vertex- 
correction to the LO diagram and the other from inserting a chain into one of the fermion 
loops of the chain in the LO diagram. Each contribution individually gives (n — 1)! for 
large n, but this leading terms cancels in the sum of both by the cancellation discussed 
at the end of Sect. 3. As a result we find 



IPR.NLO ri^co _ J_ C 



2 



const-/?^(n-2)!7^«)7'^, (5.42) 



where the constant is determined by the 1/n-correction to the asymptotic behaviour of 
the leading order vacuum polarization, see e.g. Sect. 5.1. Compared to the IPI diagrams, 
the IPR ones are suppressed by for large n. 

To compare this with the solution to the renormalization group equation, we expand 
eqs. (|5.35|) and (|5.36| ) in 1/iV/ and note that up to 1/n-corrections, the large-order 



behaviour of the four-point function is given by 



4 

^ {Ei{a) 7^75 J^J5 + E2{a) 7^ 7^^) . (5.43) 



With 



1 99 99 , 



we get 

Si(a) = ^/3Jn!a"+^^i ^2(0) = ^ /^o (™ - 2)!a"+^ ^2 , (5.45) 

n n 

where the 1 / n^-suppression of E2 follows from 71 = —2 at leading order. By comparison 
with the NLO calculation one determines 

E, = -^e^, (5.46) 
2Nf ' ^ ^ 

while E2 is related to the 1/n-correction to the leading order vacuum polarization. 

At NNLO we separate again the IPI from the IPR diagrams. All terms that involve 
7/i®7^ have been evaluated in Sect. 4.1.2. All other IPI diagrams that did not contribute 
to the vertex function in Sect. 4.1.2 are proportional to 7^75(8)7'^ 75 in large orders and do 
not yield logarithmic enhancement. We can summarize their contribution by a constant 
Dbox in order to obtain 



,lPI,NNLO 



7^1 n 



^(1)) + 



Nf 



7m75 ® 7''75 
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+ 



2Ni 



(lnn-V(l)) + ^|7M®7'' 



(5.47) 



with Cbox as in Sect. 4.1.2. The IPR diagrams involve insertions of the next-to-leading 
order results for the vertex function and vacuum polarization collected in Sect. 4. There 
is only a partial cancellation between vertex and vacuum polarization insertions and the 
IPR diagrams are not suppressed by l/ri^ in comparison to the IPI ones any more: 



„lPR,NNLO "^oo 



97r 



Nf 



^ (Inn-V'(l)) 
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(5.48) 



In the sum of IPR and IPI contributions only the last term survives in the 7^ ® 7^^- 
structure. Expanding eqs. ( |5.35[ ) and ( p.36| ) to NNLO, we have 



El (a) = 5]/3o"n!a"+i 



El ( 1 + ^Inn 
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11 

'Wf- 
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(5.49) 



neglecting 1/n-corrections. Given Ei above, E2{a) is completely determined and co- 
incides with the 7^ ® 7'^-term in the explicit calculation, the sum of eqs. ( p. 47] ) and 
(|5.48| ) - a highly non-trivial consistency check. The comparison of the 7^75 ® 7'^75-term 
determines the l/A^|-corrections to Ei. 

Although the initial conditions can only be determined order by order in 1/Nf, the 
n-dependence in large-orders is determined already by the anomalous dimension matrix. 
Since 72 > 71, the final result for four-fermion scattering can be written as 



X 2 



97r 



l + Ol 



11 



7/^75 ® 7'^75 



4N 



f 



(5.50) 



For Nf = 1, (3i/ (3q + 72 = 4.48.... In the above equations all terms {\n.n/Nf)^ are 
summed to all orders in 1/A^/ and corrections are of order 1/n and 1/Nf to the normal- 
ization as indicated. 

Vacuum polarization. The asymptotic behaviour up to l/n-corrections follows 
from insertion of O^^ into the photon two-point function, which results 



ts m 



[-^) {Qh^^u - QM ^4 (a) 



(5.51) 



Again we check that the renormalization group equation sums correctly the Inn-terms 
in the finite-order 1/iVj-result. Since we did not compute 1/n-corrections to the leading 
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asymptotic behaviour, we can drop the terms containing rp^ ~ (a in eq. (|5.38| )) 
in comparison to those with rC'^. Furthermore, since A2 ~ V^/ -^1 ^'^d -^2 are of the 
same order in 1/Nf, we may also neglect all terms that contain A2-E2 in the remaining 
terms. With these simplifications 



n ^ 

/ (2Nf + 1^ 2 



1/2 



11 



+ 11 7^3 73iAi- 732ir 



6 J 72(1 + 72) 
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(5.52) 



Notice the term n'^^/72, which, upon expansion in 1/Nf, produces 

n — independent + In n + O ( — — ) (5.53) 

with 77,-independent terms proportional to Nj. These terms would be absent, had we 
expanded the renormalization group equation first in and then integrated it. In 

this case we would have obtained Inn for rO'^ 1^2-, when it occurs in £^3(0;) and Inn + 1, 
when it occurs in Eijyoi). (The additional comes from expansion of 1/(1 + 72) in 
£'4(0), when the expansion in is performed after integration of the renormalization 
group equation.) The terms proportional to positive powers of N j are indeed irrelevant 
and can be absorbed into a redefinition of i?3. In the following this redefinition will be 
understood implicitly, so that we may substitute rf'^ /'^2 by Inn (in E^{a)) and Inn + 1 
(in E^{q)) to the order in 1/Nf, where we can compare with explicit calculation. Finally, 
we obtain 



E^ia) = Poril n ( 1 + ^ Inn 



E4 H ^nE-s + - — ^ n(lnn + 1 



(5.54) 



where the leading order result for Ei has been used. Explicit calculation, including the 
light-by-light scattering diagrams, yielded 
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(5.55) 



as coefficient of ia'^'^^. Comparison of the previous two equations shows that the Inn- 
term is indeed correctly reproduced (recall the factor of two in eq. (|5.51| )) and the others 
can be accounted for, if 

E, = -L.e^ (5.56) 
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and the terms proportional to n in square brackets, including Cbox from light-by- light 
scattering, are absorbed into E3. The correct asymptotic behaviour of the vacuum 
polarization is 



const X ^Jn! n2+^i/^o+72 ^ 



(5.57) 



the constant again being determined only as expansion in l/iV/- 

Our solution for -£^4(0) that determines the asymptotic behaviour of the vacuum po- 
larization differs from , because we consider a different theory, QED with Nj fermions 
of identical charge. In [|T^, the condition J2iQi = ^ on the charges of the fermions has 
been imposed, so that the light-by-light contributions to the vertex and vacuum polariza- 
tion vanish. This condition modifies the anomalous dimension matrix and consequently 
the solution for Ei{a) for i = 3,4. As far as we can tell, our result would coincide 
with Q, if the same condition were imposed, although our derivation and formalism is 
technically different (but physically equivalent). 



Vertex function. Since all Ei are now fixed (to leading order in 1/Nf), we obtain 
a prediction for the vertex function without free parameters. For the three-point func- 
tion, the IPR diagrams with an insertion of O4 into the 'external' photon line are not 
suppressed for large n, since E^^a) / E3{a) ~ n. The large-order behaviour is given in 
terms of E^i^a) and £'4(0;) by 



9-^ {Q^lti- 



[E^{a)-g^E^{a)]. 



(5.58) 



After expansion in 1/Nf, the terms involving E^ for i = 1, 2, 3 almost cancel in the com- 
bination -£^3(0;) — g'^E^la), the difference arising entirely from the different interpretation 
of rC^ /'')2 for E^{q) and E/^{a), when expanded in 1/Nf, as mentioned before. We are 
left with 



81n 



(5.59) 



An analogous cancellation between IPI and IPR diagrams occurs in the explicit calcu- 
lation to NLO in 1/Nf, for which we obtain 
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The second line comes from the vertex function in Sect. 4.1.3 and the third from insertion 
of the vacuum polarization (Sect. 4.2.2) into the external line. The sum is in agreement 
with the previous result. To obtain the IPR vertex function, one would have to amputate 
the external photon line and the corresponding contribution to the large-n behaviour that 
comes from the photon vacuum polarization. 

6 Conclusion and outlook 

We have described a systematic approach to factorial divergence in large orders of per- 
turbation theory generated by large-momentum regions of loop integrations. In this 
approach we reorganize the full perturbation series in terms of diagrams with a fixed 
number of chains, equivalent to an expansion in 1/Nf. We can then exploit a standard 
technique applied for analysis of UV/IR divergences and singularity structure of ana- 
lytically regularized Feynman amplitudes to classify the UV renormalon singularities of 
the Borel transform for any given class of diagrams. The problem of finding the over- 
all normalization K of the large-order behaviour reduces to calculating the residues of 
these singularities. The coefficient of the strongest singularity can be found algebraically 
by repeated contraction of one-loop subgraphs, but provides information on anomalous 
dimensions only. A non-trivial contribution to K requires calculation of part of the 
next-to-leading singularity. However this part is local so that it can be obtained with 
the help of infrared rearrangement. Still, diagrams with an arbitrary number of chains 
contribute to the overall normalization. 

Using the similarity of our UV renormalon calculus with the analysis of usual loga- 
rithmic UV divergences, we have shown that the contributions from UV renormalons to 
the large-order behaviour are naturally characterized as insertions of higher-dimensional 
operators into Green functions, a hypothesis originally forwarded by Parisi These 
insertions assume the form J2iEi{a)Goi{l^ Q^); where the Ei are independent of external 
momenta and account for factorial divergence from loop momentum regions k^q, where 
all loop momenta are much larger than the external momenta. On the other hand higher- 
order corrections in a to the Green function with operator insertion GqXi'i ^) ^ske into 
account subleading contributions in n for large n from loop momentum regions, where 
at least one of the loop momenta is of order of the external momenta. 

This systematic organization comes at a price, the introduction of an artificial expan- 
sion parameter 1/Nf. In particular, diagrams with a larger number of chains (and thus 
suppressed in 1/Nf) are not suppressed for large n, the order of perturbation theory. 
Quite to contrary, they are typically enhanced by a factor of In n per chain, so that any 
finite order in the l/A^/-expansion does not provide the correct asymptotic behaviour in 
n in the full theory for any Nf. However, just as renormalizability of the theory allows 
us to go beyond a finite-order expansion in a using renormalization group equations, the 
fact that UV renormalons can be compensated by counterterms of higher-dimensional 
operators to any order in 1/Nf allows us to transcend a fixed-order expansion in 1/Nf, 
using low-order results only. The application of renormalization group ideas then leads 
to resummation of (ln?T,/A^j)'^-terms and the correct large-n behaviour is found in terms 
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of anomalous dimensions of dimension-six operators (see eq. ( |5.57| ) for the photon vac- 
uum polarization). The overall normalization is determined by matching order by order 
in 1/Nf. The behaviour of the 1/A^j-expansion for this constant remains a problem that 
can not be addressed within this framework. Because the overall normalization remains 
undetermined for all practical purposes in the realistic situation {Nf ~ 1), a modified 
QED without UV renormalons (quite similar to the suggestion of |^) - and therefore 
potentially Borel summable - resembles in effect a non-renormalizable theory. For each 
higher-dimensional operator an unknown constant Ei has to be introduced. Of course, 
this is just the statement that for all practical purposes QED can be considered as an 
effective theory below the Landau pole, not withstanding its hypothesized 'triviality'. 

On the technical side, we have somewhat glossed over the details of two issues. The 
first is the calculation of 1/ n-corrections to the leading asymptotic behaviour. Here one 
could also use IR rearrangement, provided one properly subtracts the non-local terms 
that arise from subdivergences. This procedure can lead to a non-trivial interplay of 
analytic and dimensional regularizations. A similar interplay affects the second issue 
of how to implement the correct MS subtractions (or any other) implied to obtain the 
renormalized perturbation series in the first place. In QED this issue concerns only the 
overall subtractions for the photon vacuum polarization. While these subtractions do not 
involve factorial divergence at the level of a single chain, their insertion as counterterm 
into more complicated diagrams contributes to factorial divergence at the level of 1/n- 
corrections to the leading asymptotic behaviour of diagrams with two or more chains. 
Thus both issues must be considered when computing l/ra-corrections, but do not affect 
the status of the overall normalization K, which is independent of subtractions up to 
the trivial factor e*". 

There exist several obvious extensions of this work. The techniques developed here 
are applicable without modification to the next UV renormalons {u = 2,3, . . .), which 
are related to operators with dimension larger than six. However, one has to cope with 
a significant increase in complexity, because multiple insertions of dimension six (and 
in general, lower-dimensional) operators contribute and the Borel transform is a multi- 
valued function. When applying the formalism the Heavy Quark Effective Theory (to 
be precise, an abelian version of it), owing to the different power-counting for the static 
quark propagator, one can obtain singularities at half-integer values of u. A special 
feature of this theory is that the self-energy is linearly UV divergent. If one repeats the 
analysis of the heavy quark self-energy along the lines of Sect. 3, one finds that the 
terms iVf and Nf, not computed in [^, add to zero. 

A natural extension would include infrared renormalons. The corresponding poles 
at negative u can again be classified in terms of singularities of analytically regularized 
integrals, now from regions of small loop momentum flow through chains. The modifica- 
tions of our formalism arise only from the different notions of irreducibility of subgraphs 
in the infrared and ultraviolet. Thus IPI (UV-irreducible) subgraphs in Sect. 3.1 should 
be replaced by IR-irreducible subgraphs (UV) forest by IR forests etc. For the 
general structure of large-n behaviour associated with diagrams with multiple chains, we 
expect little difference from the ultraviolet case. In the language of counterterms in the 



Lagrangian, IR renormalons would be interpreted as non-local operators |36], since IR 



renormalons are local in momentum space (while UV renormalons are local in coordinate 
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space) . 



The ultimate goal is still a diagrammatic understanding of renormalons in non-abelian 
gauge theories, that is QCD. Given that even in QED with Nf ~ 1 there is little to say 
about the overall normalization of large-order behaviour, there is little to hope that one 
could obtain more than the exponent b in the expected behaviour 



Already this form poses the challenge to understand, on a diagrammatic level, how the 
contributions from many diagrams with many loops combine to produce the first co- 
efficient of the non-abelian /3- function, which, contrary to QED, is not related to 
loop-insertions into the gluon propagator. Although the physics described by non-abelian 
gauge theories changes discontinuously with Nf for sure, the structure of Feynman dia- 
grams does not, so that we may still try to approach the question from large Nf, provided 
we consider the expansion in Nf as formal. Writing (3^^ = /^o + where (3o denotes 

the abelian part, the previous equation is expanded as 



so that for each chain (defined as in the abelian theory as gluon propagator with fermion 
loops) we need an additional enhancement of n (rather than In n) that combines exactly 
to to all orders in In QCD, contrary to QED (where these contributions 

cancel between various diagrams due to the Ward identity), one also has to take into 
account the possibility of first inserting a (usual) UV counterterm and then an UV 
renormalon counterterm, in addition to the situation discussed at length where one first 
inserts a UV renormalon counterterm and then a (usual) UV counterterm, related to 
renormalization of a dimension-six operator. A more detailed investigation reveals that 
the first contributions combine to produce the term 5(3^^71 above, while the others are 
interpreted in terms of dimension-six operators just as in the abelian theory. 
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A Singularity structure of analytically regularized 
Feynman integrals 

A.l ^-representation 



A general theorem |]T8[ regarding the analytical structure of an analytically regularized 
Feynman integral -Frl^jIZi; A) states that it is a meromorphic function of the regulariza- 
tion parameters. 

A standard way to prove this property and to get more concrete information about 
singularities with respect to A is to apply the well-known a-representation of the Feynman 
integrals^ which is obtained by rewriting propagators as 



(p2 -)- 77j,2) 



1 



and performing GauB integration in the loop momenta in the integrand of the integral 
in ai, . . .ai- Since the parameters r; can be taken into account by a redefinition of A; 
we suppose now that r/ = for all /. 

After that the analytically regularized Feynman integral takes the form 

Fr{q,m;y) = C(A) / daY\af'D{a)-^Z{q,a)expl-A{q,a)/D{a) -Y^mfai] 

I ~ [ ~ 1=1 J 

= [ daY[af'Iiq,rn,a;X) , (A.2) 
•' I 

where 

Z{q,a)=l[Zi{d/dui) e^p{{B{q,C,a) + K{q,^,a))/D{a)}\^^^^ . (A.4) 

Here ^ = (^i, . . . , a = (ai, . . . , ul), da = dai . . . da/,. The integration is over the 
domain of non-negative parameters a, and D, A, B, K are homogeneous functions that 



are constructed for the graph F according to well-known rules (see, e.g., ^ p^j). 
The homogeneity degrees of the functions D{a) and A{q,a) with respect to the set of 
variables a are respectively h and h + 1. The function A{q^a) is quadratic with respect 
to q. 



^^For simplicity, we consider Euclidean Feynman integrals. The analytical properties in A are the 
same in Minkowski and Euclidean spaces. 
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A. 2 Sectors 



The problem of characterizing the singularities in A can be reduced to resolving singu- 
larities of integrands in the a-representation which is a problem of algebraic geometry. 
A natural approach is to locally introduce new variables in which the original compli- 
cated singularities are factorized. In renormalization theory however one applies a simple 
change of variables to achieve this goal. First, the integration domain in ( [A.2| ) is de- 
composed into subdomains that are called sectors. Second, in each sector new (sector) 
variables are introduced. In these variables the singularities of all the functions involved 
are factorized so that the problem reduces to power counting in the sector variables. 

We will not deal with IR problems so that it is sufficient to introduce sectors and 
sector variables associated with IPI subgraphs and corresponding to maximal generalized 
forests. Let us call a subgraph UV-irreducible if it is either IPI or consists of a single line 
which is not a loop line. A set F of UV-irreducible subgraphs is called a generalized forest 
if the following conditions hold: (a) for every pair of IPI subgraphs 7,7' G -F one has 
either 7 C 7', or 7' C 7, or 7 and 7' are disjoint (with respect to the set of vertices); (b) if 
7i, . . . , 7j G F are pairwise disjoint with respect to the set of fines then the subgraph Uj7i 
is one-particle-reducible (IPR). Generalized forests which do not involve single lines are 
nothing but forests in the common sense. Thus every generalized forest F is represented 
as the union Fh U of a forest Fh consisting of IPI elements and a family Fr of IPR 
single lines. 

Let JF be a maximal generalized forest so that for any 7 which does not belong to 
JF the set U {7} is no longer a generalized forest Let a be the mapping defined by 
(7(7) G 7 and (7(7) ^ 7' for any 7' C 7, 7' G JF so that (7(7) is the fine that belongs to 7 
and does not belong to subgraphs smaller than 7. 

Let us define the following family of domains (sectors) [l^, 
the maximal generalized forests of P: 

= < "(t{7) V/ G 7 G J'^} . (A-5) 

The intersection of any two distinct (corresponding to different generalized forests) sec- 
tors has zero measure and the union of all the sets Ujf "^^j^ is the whole integration domain 
in the a-representation. 

Let F^{q,rn;X) be the contribution of a sector Vjr to the given Feynman integral. 
Let us introduce the following (sector) variables: 

ai= n ^- (A-6) 

The corresponding Jacobian equals Yl^,zjrt^^"'^~^ . The inverse formulae are 

t-y = acr(7)/"CT(7+) ) (A. 7) 

where 7+ denotes the minimal element of JF that contains 7 (we put ao-(r+) = !)• 

Due to factorization properties of the homogeneous functions of the a-representation, 
Fr can be represented as 



18, EOI associated with 
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Fr^(g,in;A)= ri?(^)-['^™-Mir R /^^^^^-I^^^^/^l-MVl^,!",*; A) , (A.8) 

where / is an infinitely differentiable function, the square brackets denote the integer 
part of a number, and 

m = j2xi. (A.9) 

ley 

Thus the UV convergence analysis and therefore the analysis of the analytic properties 
in A is characterized with the help of properties of one- dimensional distribution 
which is a meromorphic function with respect to A, with simple poles at the points 
A = — 1 — n, n = 0, 1, 2, . . . where 

4 = i^^<5W(x) + 0(l). (A.10) 



A. 3 Taylor operators 

To analyze analytical structure of the Feynman integral with respect to complex param- 
eters A let us introduce standard Taylor operators and present them in various forms. 

Let be the operator that picks up the terms up to the A^-th order of the Taylor 
expansion in the corresponding set of variables. Then the formal Taylor expansion in 
masses and (independent) external momenta is described in the a-parametric language 
as 



K=l 



= lda^af'T,^K^'^(^^+^(^h{q,m,a'-,X)l^^, (A.ll) 

where a'l = H?ai. 

Note that this is just a formal expansion operator in the sense that it would certainly 
generate IR divergences when naively applied to the Feynman integral. But we shall 
usually apply such operators not to the whole integrals but to its sector contributions. 
These operators are by definition applied to the integrand of the a-integrals involved. 

Let 7 be a subgraph of F. Let us perform its mass and momentum expansion by 
the operator and then insert the result into the reduced graph F/7. We denote the 
insertion of a polynomial V in external momenta of the subgraph 7 into the reduced 
diagram by Fr/-y o V. Thus we are dealing with Fr/^ o M-yFr^ which is by definition the 
action of the operator My associated with 7 on the whole Feynman integral Fp. 

To describe this procedure in the a-parametric language let us remember that the 
tt-representation is obtained by performing integrals over the loop momenta. If one first 
integrates over thee loop momenta of the subgraph 7, one obtains the a-representation 
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for the Feynman integral associated with the subgraph. After using formula ( A.ll|) and 
continuing integration over the rest loop momenta (in fact associated with the reduced 
graph r/7) one arrives at the relation 



M^^Fr(£,m;A) = (Fr/7 o A/^F^)(g, m; A) 



,rn,a';X) 



K = l 



(A.12) 



where a'l = K^ai if 7 and = a; if / G 7. 

These simple arguments would be rigorous if one could control convergence of the 
integrals in the loop momenta and in the a-parameters involved and easily prove that one 
may change the order of the integration. It happens that it is simpler to provide a direct 



proof of ( |A.12| ) — see [20|| (Lemma 5). More details can be found in |^| (Lemma 11.1). 



Let 7 G JF. Note that by ( [A. 6] ) multiplication of all ai with / G 7 by is equivalent 
to the replacement n'^t^. Therefore the action of the operator on the sector 

contribution takes the form 



(A.13) 

where / is the function that enters (|A.8|) . 



A. 4 Analytic structure in A 

Representation ( |A.^ ) shows that the sector contribution {q^uv^X) to the analytically 
regularized Feynman integral i^r(9,2Zi;A) can be represented as a 

where the function gjr is analytical in a vicinity of the point A = 0. The whole Feynman 
integral is therefore looks like the sum ( |2.20|) of the terms (|A.14|) over all the maximal 



forests of F. 

The leading singularity of the given Feynman integral is by definition the sum of terms 
( [A. 141 ) with the maximal number of the factors A (7) in the denominator. According to 



([A.14| ) this number is equal to or less than the maximal number of divergent subgraphs 



that can belong to the same forest. For simplicity, let us now consider only the vicinity 
of the point A/ = 0. 

To calculate the leading singularity let us observe that the singular factor 
originates from the integral over the sector variable in ( |A.8| ) . Due to ( |A.1CI| ) , evaluation 
of the residue at this pole amounts to the action of the Taylor operator T^^'^'^ which in 
turn is equivalent to the action of the operator M/^, with N = a;(7), given by the 



right-hand side of ( |A.12| ). Using this relation we observe that the problem reduces to 
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evaluation of the residue of the Feynman integral -^^(A) with respect to A (7) and insertion 
of the result into the reduced graph. 

Consider, for example, a situation when a given maximal forest contains only two 
divergent subgraphs 7 and F. Let Kx be the operator that picks up the pole part of 
Laurent series in A. Then we see that the leading singularity LS{Fy) is given by 

KxiT) (i^r77'(Ar/7,{A(7)}) oi',(,)Ff ^)(A)) , (A.15) 



Remember that A(7) is given by (|A.9|) . Here jF/7 is a generalized forest of the reduced 



graph (it is obtained from by 'dividing' each element 7' with 7 C 7' by 7) and ^^(7) 
is the generalized forest of 7 which is the 'projection' of JF to 7. (Note that there is one- 
to-one correspondence between generalized forests JF that satisfy 7 G JF and generalized 
forests of 7 and F/7.) By F^j^i- ■ ■ , {A}) we denote the given sector contribution regular- 
ized by inserting the factor cta{j+) = '''7+^f7+)+ • • • i^^o ^^e integrand. Finally, F'^^''\x) 
differs from -F!^'-''''' ( A) by the reduction of the region of integration to ai < a;o-(7+), I G 7. 
Let us now observe that 

i^A(7)^f "^(A) = k^.^F^^^HX) = k,F^^^\X)\^^^^^^ = _l_resAF^^W(A) . (A.16) 

where it is implied that F:^^"'\X) is regularized by introducing the regularization param- 
eter A into any line of 7. Similarly, we have 

i^A(r) {F^/;arH, {A(7)}) oV) = K, (f,7;(A) o v) ^^^^^^ ^ ^res. (f,7;(A) o v) , 

(A.17) 



Here V is the polynomial (|A.16| ) that is inserted into the reduced graph. 



Let us represent a maximal forest JF as JF/^ U JF^ = J-'h,div U ^h,conv U JF^ where the 
subscripts 'div' and 'conv' denote IPI elements respectively with uj{-y) > and ^(7) < 0. 
Let maxj^ — '^o be the maximal possible number of divergent subgraphs that can 

belong to the same maximal forest.]^ 

Using straightforward generalization of the arguments used in the above example with 
two divergent subgraphs and performing summation over all the maximal forests of the 
given graph we come to the following simple representation for the leading singularity: 

LS{Fr)= n 3^(resAF,/^_(A)) . (A.18) 

Here 7_ is the set of maximal elements 7' G JF with 7' C 7. Each factor TesxF^/^_ (A) is a 
polynomial with respect to external momenta and internal masses of 7/7.. It is implied 

-'^^Note that a single line with the propaga tor l/(p^)^ is formally UV-divergent. However the corre- 
sponding pole is cancelled by a F- function in (|A.2| ) . Therefore it is sufficient to consider the usual forests 
rather than generalized forests when calculating the number of singular factors that enter the leading 
singularity. 
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that these factors are partially ordered and before calculation the residue TesxFj/.y_{X) 
all polynomials associated with the set 7_ are inserted into this 'next' reduced diagram. 
Remember that when calculating these residues it does not matter into which line the 
regularization parameter A is introduced. 
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